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PREFACE 


The  research  presented  in  this  paper  represents  the  bulk 
of  the  work  performed  on  the  numerical  integration  of  simple 
non-linear  atmospheric  models  by  the  Numerical  Prediction  Project, 
a  group  jointly  supported  by  the  Geophysics  Research  Directorate 
and  the  Air  Weather  Service,  during  the  period  from  the  project's 
activation  in  February,  1953  to  August,  195^.  Members  of  the 
project  are  connected  with  the  Directorate's  Atmospheric  Analysis 
Laboratory . 

Partly  to  determine  the  operational  feasibility  of  the  models 
then  available  for  numerical  forecasting  pr-noses  and  partly  to 
lay  the  foundation  for  a  systematic  researc..  rogram  in  numerical 
prediction,  it  was  decided  to  examine  the  performance  of  the  fam¬ 
iliar  barotropic  model  and  of  a  simple  baroclinic  model  under  a 
wide  variety  of  synoptic  conditions.  The  month  of  January,  1953* 
was  selected  for  this  test,  comprising  a  continuous  series  of  60 
cases . 


The  efforts  of  the  members  of  the  project  were  then  focused 
on  the  development  of  a  suitable  baroclinic  model,  the  plotting, 
analysis  and  tabulation  of  the  required  synoptic  data,  the  design 
and  construction  of  the  programs  and  codes  for  solution  on  a  high¬ 
speed  computer,  and  finally  on  the  analysis  and  interpretation  of 
the  results.  The  overall  organisation  of  this  effort  was  directed 
by  Major  P.  D.  Thompson,  who  was  also  responsible  for  the  develop¬ 
ment  of  the  thermotropic  equations  tested.  The  analysis  and  tab¬ 
ulation  of  the  synoptic  charts  was  supervised  by  Mr.  C.  F.  Jenkins, 
while  the  calculation  and  summary  of  the  statistical  properties  of 
the  forecasts  was  supervised  by  Dr.  I.,  3.  Pocinkl.  Dr.  W.  L.  Gates 
assisted  Major  Thompson  in  supervision  or  preparation  of  the  numer¬ 
ical  forecasts,  and  was  responsible  for  development  of  the 
thermotropic  code.  The  code  for  solution  of  the  barotropic  model 
was  written  by  Major  H.  A.  Partner.  The  successful  completion  of 
the  tests  was  materially  aided  by  other  members  of  the  project 
during  this  period,  including  Mr.  L.  Derkofsky,  Mr.  E.  J.  Aubert, 
Major  J.  F.  Blackburn,  Lt .  M.  E.  Stern,  Mr.  J.  J.  Pazniokas, 

Lt .  Patricia  Hayes,  and  Sgts.  D.  Casey  and  D.  Grennan.  Special 
recognition  is  due  Mr.  E.  A.  Bertoni  for  his  preparation  of  the 
figures  in  this  report.  Mr.  W.  S.  hering  and  Mr.  W.  D.  Mount  of 
the  Atmospheric  Analysis  Laboratory  assisted  in  verification  and 
synoptic  summary  of  the  forecasts. 
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PREFACE  (Continued) 


The  details  of  the  theoretical  development  and  a  state¬ 
ment  of  the  overall  results  have  been  presented  in  a  separate 
paper.  This  paper  has  been  prepared  to  present  the  detailed 
analysis  of  the  results  of  this  series  of  numerical  forecasts, 
in  the  hope  that  they  may  serve  as  a  standard  of  comparison 
with  other  methods  of  prediction  and  allow  the  improvement  of 
such  techniques  of  dynamical  prediction  through  the  focusing 
of  further  research  on  their  Inadequacies.  In  particular,  it 
is  hoped  that  the  statistical  and  synoptic  analyses  o'  in. 
results  will  permit  an  evaluation  of  these  techniques  by 
synoptic  meteorologists  unfamiliar  with  the  theoretical  aspects 
of  numerical  prediction.  To  further  these  purposes  the  actual 
24 -hour  numerical  forecast  charts  for  1500  G.C.T.  of  each  day 
of  the  series  are  presented  in  the  appendix  to  this  paper. 
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ABSTRACT 


Following  a  resume  of  the  theory  of  thermotropic  flow, 
a  simple  baroclinic  model  in  which  the  direction  of  the 
thermal  wind  is  assumed  invariant  with  height,  a  discussion 
of  the  methods  employed  for  the  numerical  integration  of 
this  model  and  the  barotropic  model  is  presented. 

On  a  finite-difference  grid  of  4l4  points  covering 
the  United  States  and  immediately  surrounding  regions,  a 
series  of  sixty  comparative  24-hour  forecasts  during  Jan¬ 
uary,  1955,  at  the  50G  and  1000  mb  levels  was  obtained  by 
relaxation  methods.  The  median  correlation  coefficients 
between  the  forecast  and  observed  24-hour  height  changes 
were  0.6j  for  both  thermotropic  and  barotropic  models  at 
500  mb  and  0.69  for  the  thermotropic  model  at  1000  mb.  In 
comparison  with  a  method  of  pure  interpolation,  the  numer¬ 
ical  forecasts  are  shown  to  display  a  positive  "skill" 
toward  the  center  of  the  forecast  region  where  the  influence 
of  the  lateral  boundaries  is  smallest.  By  normalising  the 
root -mean-square  forecast  errors  to  allow  for  the  normal 
latitudinal  variation,  the  Rocky  Mountains  are  found  to 
exert  a  marked  influence  on  the  forecasts  at  both  500  and 
1000  mbs  over  the  south  central  United  States. 

From  a  synoptic  point  of  view,  the  numerical  forecasts 
are  found  to  compare  favorably  with  conventional  forecasts 
for  the  same  period,  although  they  appear  to  introduce  a 
small  but  systematic  tendency  to  move  fully  developed  dis¬ 
turbances  too  slowly.  This  error  is  felt  to  stem  from  the 
truncation  errors  of  the  finite-difference  schemes  employed. 

Recommendations  for  further  research  to  reduce  the 
several  sources  of  error  and  to  extend  the  physical  basis 
of  the  model  are  made. 
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RESULTS  OF  NUMERICAL  FORECASTING  WITH  THE  BAROTROPIC 


AND  THERMOTROPIC  ATMOSPHERIC  MODELS 


1.  General  Introduction 

During  1952  a  number  of  workers  succeeded  in 
devising  several  alternative  baroclinic  models  (approxi¬ 
mate  formulations  of  the  general  theory  of  baroclinic 
flow)  that  are  workable  from  the  computational  standpoint, 
and  which  contain  some  of  the  ingredients  essential  to  the 
intensifications  of  weather  disturbances.  These  models 
include  the  two-parameter  model  of  Eliassen  ,2  the  so-called 
"2  l/2-dimensional"  modeLof  Eady,^  the  equivalent-baro- 
clinlc  model  of  Thompson,  the  two-level  model  of  Charney  , 
and  Phillips, 5  and  the  two-layer  model  of  Sawyer  and  Bushby. 
All  are  substantially  equivalent. 

These  methods  had  been  applied  in  a  few  selected 
cases  of  rapid  development  with  very  encouraging  results; 
the  coefficients  of  correlation  between  the  observed  and 
predicted  behavior  of  the  height  of  the  500  mb  surface 
averaged  around  0.7,  and  indicated  that  numerical  methods 
might  be  successfully  applied  in  situations  of  strong 
cyclogenesis.  Since  these  results  could  be  assumed  to  hold 
in  less  extreme  weather  situations,  it  was  generally  agreed 
that  the  methods  of  numerical  weather  prediction  showed 
promise  of  producing  a  significant  increase  in  the  accuracy 
of  short-range  forecasts.  These  methods  might  also  be  ap¬ 
plied  on  an  operational  basis. 

I*,  must  be  emphasized,  however,  that  even  the 
simplest  of  these  baroclinic  models  had  been  applied  only 
in  a  few  specially  selected  and  rather  unusual  cases.  It 
t none  fore  appeared  reasonable  that  several  methods  of 
numerical  prediction,  varying  as  to  generality  and  com¬ 
plexity,  should  be  subjected  to  extensive  comparative  tests. 

The  present  series  of  tests  was  consequently 
organized  with  the  following  general  objectives  in  view: 
il)  to  establish  the  over-all  accuracy  of  several  existing 
methods,  by  applying  them  to  a  standard,  representative 
sample  of  weather  situations,  and  therby  lay  the  foundation 
'‘or  a  continuing  program  of  research  in  numerical  fore¬ 
casting,  and  (2)  to  develop  and  standardize  procedures  for 
producing  forecasts  on  a  regular  basis  from  existing  or 
anticipated  methods  of  numerical  prediction. 
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This  research  program  involved  the  performance  of 
a  comparative  series  of  numerical  forecasts,  based  on  differ¬ 
ent  methods  of  prediction  but  on  the  same  observational  data. 

It  was  planned  to  compute  12-  and  24-hour  forecasts  from  60 
consecutive  sets  of  initial  data  spaced  12  hours  apart,  in 
the  expectation  that  errors  could  be  more  easily  isolated  and 
analyzed  in  a  consecutive  series  than  in  a  series  chosen  at 
random  from  a  larger  sample.  The  month  selected  for  study 
was  January  1953 ,  when  radiosonde  and  rawinsonde  observations 
over  the  United  States  were  quite  dense,  and  when  a  fairly  rep¬ 
resentative  variety  of  weather  conditions  prevailed  over  the 
continental  United  States  and  environs. 

After  preliminary  study,  it  was  decided  to  compute 
forecasts  from  the  equations  for  the  barotropic  model  and  for 
the  so-  ailed  "thermotropic"  model,  a  vertically-integrated 
two-pa. -meter  model  summarized  in  Section  2.  Comparisons  be¬ 
tween  these  two  parallel  series  of  numerical  predictions,  the 
changes  actually  observed  to  occur,  and  the  forecasts  prepared 
by  conventional  methods  would,  of  course,  comprise  an  essential 
first  step  toward  the  objectives  of  the  program.  These  com¬ 
parisons  from  both  a  statistical  and  a  synoptic  viewpoint  are 
presented  in  Sections  3  and  4,  wherein  the  present  technique 
and  theory  of  numerical  prediction  is  evaluated  both  as  a  re¬ 
search  tool  and  as  a  forecasting  technique. 


2 .  The  Numerical  Integration  of  the  Thermotropic  Model 

2.1  Resume  of  the  Thermotropic  Equations 

The  derivation  of  the  thermotropic  qua t ions  for  an 
adiabatic,  frictionless  atmosphere  in  hydrostatic  and  quasi - 
geostrophic  balance  has  been  described  in  detail  by  Thompson 
and  Gates, *  but  will  be  summarized  here  in  order  to  provide  an 
introduction  as  well  as  a  background  to  the  detailed  discussion 
of  results  in  subsequent  chapters.  This  so-called  thermotropic 
model  is  a  baroclinic  two -parameter  model  for  vertically-inte¬ 
grated  flow,  of  which  the  basic  meteorological  assumption  is 
that  of  an  isogonal  thermal  wind,  i.e.  the  direction  of  the 
wind  shear  does  not  change  with  height  although  the  magnitude 
of  the  shear  may  change.  In  terms  of  the  temperature,  T,  the 
thermotropic  or  two -parameter  model  may  be  charac cerized  by 
the  equation, 


vr  =  F(f>)  vf 
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where  p  is  the  isobaric  gradient  vector,  and  where  T 

-  bo  "V7:^  ,  with  the  function  F(p)  expressing  the 

(magnitude)  variation  of  the  thermal  wind  vector.  This 
characterization  of  the  atmosphere's  baroclinity  is  es¬ 
sentially  that  embodied  in  the  two -parameter  or  2  1/2- 
dimensional  models  of  Eady,  Sawyer  and  Bushby,  and 
Eliassen;  it  may  also  be  shown  to  be  the  essential  basis 
of  the  two-layer  baroclinic  model  of  Charney  and  Phillips. 

Neglecting  the  vertical  advection  of  vorticity, 
by  assuming  that  the  relative  (geostrophic)  vorticity  is 
much  less  than  the  Coriolis  parameter,  and  neglecting  the 
terms  in  the  vorticity  equation  which  represent  a  "twisting" 
or  re -orientation  of  the  vortex  tubes,  the  equations  of  the 
thermotropic  model  may  be  written: 


+LT(*,Vli) 


The  bar-operator  (  J  —  fa  J%  (  )  ^  t  and  **/  Jl 
is  a  stream  function  for  isobaric  contour  height  t  ;  &  *  ZT/j. 
is  a  stream  function  for  temperature,  and  J  is  the  elevation 
of  the  terrain.  In  these  equations  >is  the  Coriolis  parameter, 
*  is  the  gas  constant  for  dry  air,  T#  is  a  representative 
temperature  at  the  isobaric  surface  near  the  ground,  f  is 
the  gravitational  acceleration,  i  is  the  northward 

rate  of  change  of  the  Coriolis  parameter  .  These  equations 
have  been  written  in  the  familiar  coordinate  system, 

with  the  positive  %  and  y  axes  directed  eastward  and  northward, 
respectively.  All  differentiations  indicated  by  the  Jacobian 
operator  /ht)(H/ty ancl  by  the  Laplacian 
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operator  V  &  =  2  +  2>  4,/dy  ,  are  to  be  carried 

out  with  the  pressure  held  fixed,  as  are  differentiations 
with  respect  to  the  time  £  . 

Equations  (1)  and  (2)  are  the  basic  dynamical 
equations  governing  the  behavior  of  the  thermotropic  model. 

The  first  three  terms  of  Eq.  (l)  will  be  recognized  as  simply 
those  of  the  barotropic  model  representing  the  local  vorticity 
changes  due  to  the  individual  conservation  of  absolute  vorti¬ 
city.  The  fourth  term  of  Eq.  (l)  represents  the  rate  of  pro¬ 
duction  or  destruction  of  mean  vorticity  by  the  advection  of 
"thermal"  vorticity  by  the  "thermal"  wind,  such  that  absolute 
vorticity  will  be  increased  in  a  disturbance  when  the  thermal 
trough  is  "lagging"  behind  the  contour  trough  -  an  easily 
recognized  and  ’well-known  synoptic  condition.  The  last  term 
in  Eq.  (l)  represents  the  mean  vorticity  generation  due  to 
the  effects  of  terrain-enforced  vertical  motions,  such  that  a 
vorticity  increase  or  cyclogenesis  tends  to  occur  in  strong 
flow  on  the  lee  side  of  elevated  terrain.  The  several  terms 
of  Eq.  (2)  can  be  similarly  interpreted  as  due  to  the  effects 
of  either  thermal  or  contour  vorticity  advections  or  to  the 
effects  of  terrain-induced  vertical  motions;  the  sixth  term 
of  Eq.  (2),  however,  represents  the  effects  of  mean  vertical 
motions  on  the  changes  of  thermal  vorticity,  and  thereby 
explicitly  introduces  the  thermodynamical  relation  for 
adiabatic  flow  into  the  system. 

Ey  introducing  the  concept  of  an  equivalence  level 
in  the  middle  troposphere,  at  which  the  mean  wind  is  equal  to 
the  vertically -integrated  wind,  the  equations  of  the  thermo¬ 
tropic  model  may  be  written  in  terms  of  the  height  of  a  selec¬ 
ted  isobaric  surface,  say  500  mb.,  and  the  height  of  an 
isobaric  surface  near  he  ground,  say  1000  mb.  With  this 
identification,  Eqs.  (l)  and  (2)  become  for  flow  over  perfectly 
flat  terrain  (  P  ) : 


»V>3( 
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+  13) 

+  tU 

-  =  0  (T) 

%  ' 


where  -fr  is  thf  500  mb  height,  and  £  is  the  500-1000  mb 
thickness.  Ire  factor  C  embodies  the  effect  of  the  height 
of  the  level  of  equivalence  in  this  model  and  is  given  by, 

ta  (l-f/f*)"'  , 

-  -  * 

where  is  the  mean  pressure,  s  -it  /a.  ,  and  f  is  the 

pressure  at  the  equivalence  level,  taken  to  be  600  mb. 


The  thermotropic  coefficients  tC  and  L  which 
appear  in  Eqs .  (1-4)  are  defined  by  the  theory  to  be  func¬ 
tions  of  pressure  alone,  and  are  given  by  Thompson  and 
Gates, i 


K  =  *  (  f  frw)  ~  t  , 


where  f(f)  is  the  thermotropic  parameter  describing  the 
vertical  variation  of  the  wind  shear  as  introduced  earlier. 
The  quantity  is  a  measure  of  static  stability  and  is 

given  by. 


- 


2  rtf  f) _ 

(fp- }>)(&,  -j.KO*x» 


with  X-  R/t )»  »  Cto  bring  the  specific  heat  of  air  at 

constant  pressure .  r 

The  simplified  dynamical  formulas,  Eqs.  (3)  and 
(4),  are  now  to  ’oe  solved  by  an  iterative  procedure  wherein 
the  system  is  regarded  as  linear  with  the  time  derivatives 
of  and  as  the  dependent  variables,  and  all  other 

• erras  regarded  a3  known  non -homogeneous  members  of  the 
equations.  Starting  from  the  known  initial  distributions  of 
500  mb  height  h  and  500-1000  mb  thickness  4.  ,  all  terms 
not  involving  time  derivatives  in  Eqs.  (3)  and  (4)  may  be 
evaluated,  a  solution  carried  out  to  obtain  0*  fit  and 
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from  the  resulting  Poisson  and  Helmholtz  equations,  respec¬ 
tively,  and  finally  a  time-extrapolation  made  from  3fc/H 
and  34/3 1  to  regenerate  £  and  4  a  short  time  later 
than  the  initial  data.  This  iterative  process  may  then  be 
repeated  cyclically  until  the  aggregate  of  the  time  extra¬ 
polations  is  equal  to  the  desired  period  of  the  forecast. 

The  actual  solutions  of  Eqs .  (3)  and  (4)  were 
carried  out  on  a  rectangular  finite -difference  grid,  and 
for  variables  expressed  in  non-dimensional  measure.  In 
the  numerical  solutions,  the  selection  «fc*s6o0^was  made, 
resulting  in  cs  0.86  when  ^  =  506  hi  jb  .  'The  thermotropic 
coefficients  K  and  L  in  these  equations  were  assigned  the 
constant  values  0.28  and  -0.l8,  respectively,*  which  were 
determined  by  fitting  the  o.bserved  winds  over  North  America 
to  an  idealized  wind  profile  for  a  series  of  ten  cases 
selected  from  January,  1953.  In  the  expression  for 
as  p;iven  by  Eq.  (8),  the  factors  in  parentheses  ..vre  found  to 
be  nearly  constant  for  the  same  series  of  cases,  with  an 
average  value  of  +  5.0  y  10“**  tt*Tz  ,  which  was  accordingly 
used  throughout  the  series  of  integrations. 

Accounting  these  modifications,  Eqs.  (3)  and  (4) 
may  be  written  in  non-dimensional  form,  with  the  omission 
of  primes  understood,  and  for  a  rectangular  grid,  as  follows: 


•  0,B6  Afj  Z>j 

-  o.et  Aqt, )(*,&)  -eJeA'aZ 

*  a, k  *  A*  a, k  *  t  u  r.'  (t,  L ) 


*For  a  linear  wind  profile,  3uch  as  that  ir.  tr 

equivalent -baroclinic  model, ^  the  values  corresponding  to 
the  thermotropic  coefficients  K  and  l,  are  1/3  and  zero, 
respectively . 


I  ^  | 

.0  coefficients  A<*  >A A  jj  and  B  •  >  are  point-wise 
vcii-iaoie  and  given  by,  'J  J  -I  J 


A'.  = 

'j 


i6*h,zi\n-  ' 


CaS 

SI 

Sih 

7L 


g1..  s 

•i 


jt.n. 


Here  H  is  the  reference  unit  of  contour  height  (such  that 
k'«  fcH"1  ,  where  is  the  non-dimensional  variable),  JL  is 
the  reference  unit  of  non-dimensional  frequency,  Jt  is  the 
mesh  constant  of  the  finite -difference  grid,  ah,  and  are 
the  mapping  or  scale  factor  and  the  cone  constant,  .  pective- 
ly,  of  the  projection  on  which  the  calculations  are  performed, 
and  ^  is  the  longitude  deviation  of  a  grid  point  from  the 
meridian  about  which  the  grid  is  symmetrical  (  <f  >  O  to  east). 
In  the  present  case.  Aft  2*2.  tv*  on  the  WBAN-1  synoptic  chart,  a 
Lambert  conic  conformal  projection  whose  scale  factor  is 

( ft /^t) ca* d L+i*  (tf/a-#)/*]*  where  &/»*  i  is  the 
ratio  of  the  pole-to-equator  distance  of  the  map  to  the  earth's 
radius,  9  is  the  geographical  latitude,  and  a,  s  6.7156  the 
cone  constant.  The  unit  H  was  taken  for  convenience  to  be 
100  ft,  SL  as  and  the  reference  unit  of  length 

was  taken  as  the  distance  between  adjacent  grid  points. 
Thus,  while  the  "physical”  coefficients  IC^L  and  of  the  ther¬ 
motropic  equations  have  been  taken  as  constants  in  the  present 
solutions,  the  "geographical"  coefficients  2,  9  and  have 
been  allowed  their  full  Of,  ^  variability. 
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The  finite-difference  operators  used  in  Eqs .  (9) 
and  (10)  are  given  by, 

A,0,v=  Ui+Ij  , 

A,  (  );j  =  l ) (j-+,  -  l  )y -I  , 

A  (  t'j  =  ( )g+i+  l\'.» j  +  (  Vj'-i  +  (  '.Vjj.  -  4  c  , 


X3U»  =  A,b6  , 

where  i  and  j  denote  measure  along  the  now  rectangular  ^  m 
axes  of  the  grid,  respectively. 

In  addition  to  the  numerical  forecasts  for  the 
500  mb  height  and  1000  mb  height  (or  500  -  1000  mb  thick¬ 
ness),  forecasts  of  the  average  vertical  motion  were  made 
from  the  first  law  of  thermodynamics  for  adiabatic  flow  in 
the  form, 


$)  4^  &>  a  Ot 


U  ~  '  T  '  '  \ p  / 

where  m  (dp /At)  ££  Vr  with  {  the  density  and 

•  When  written  in  finite-difference  form  for 
the  non-dimensional  variables  9  and  ,  the  vertical  motion 
equation  may  be  written, 


v  -  n(H  -  I,jM)}  (  uo 

where  .  The  coefficient  N  was 

found  to  be  approximately  constant,  and  was  assigned  the 
average  value-11.2  (dimensionless)  throughout  the  calcula- 
t  * ons . 

2.2  Method  of  Numerical  Integration 

With  selected  boundary  conditions,  the  method  of 
solution  consisted  of  the  determination  of  the  tendencies 
and  H/n  by  a  relaxation  method,  followed  by 
an  extrapolation  over  a  short  time  interval  tc  generate  new 
initial  data,  after  which  the  entire  process  was  repeated 
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until  forecasts  of  the  desired  length  were  obtained .  The 
finite -difference  expressions ,  Eqs .  (9)  and  (10),  are  to 
be  regarded  as  applying  to  the  finite -difference  grid  shown 
in  Fig.  4.1,  Section  4.  This  is  a^ square  grid  of  18  x  25 
points,  spaced  approximately  2  l/2c  latitude  apart  and 
covering  the  continental  United  States,  southern  Canada, 
and  the  immediately  surrounding,  territory.  This  grid  pro¬ 
bably  has  sufficient  resolution  to  faithfully  display  the 
larg  -scale  disturbances,  end  does  not  appear  to  emphasize 
small  irregularities  in  the  analysis  of  initial  conditions. 
With  the  presently  employed  method  of  determining  vorticity, 
the  complete  non-homogeneous  terms  of  Eqs.  (g)  and  (10)  are 
determined  only  on  and  within  the  third  interior  grid  points, 
yielding  solutions  for  5%  /$■£  and  /  9b  on  an  Interior 
grid  of  12  x  17  points. 


The  first  step  in  the  solution  of  Eqs.  (9)  and 
(10)  at  each  net  point  consisted  of  the  calculation  of  a 
set  of  residuals  given  by. 


07) 
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for  the  %  -  and  ^  -  equations,  respectively ,  .here  (ffc/Xi) 
and  are  a  set  of  initial  solution  estimates.  For 

*«©  ,  these  tendencies  were  taken  as  the  actually  observed 

24-hour  changes,  centered  on  the  initial  instant:  for  later 
times,  ^>0  ,  they  ..ere  taken  as  the  most  recent  solutions, 
i.e.,  the  solutions  for  the  immediately  preceding  time  stage. 
The  residues  generated  in  this  manner  .ere  probably  somewhat 
larger  for  the  first  step  lr.  the  iteration  process  than  for 
the  subsequent  ones,  with  a  consequent  increase  in  the  re¬ 
quired  amount  of  Initial  relaxation. 


The  relaxation  itself  proceeded  by  means  qf  a 
modified  form  of  the  extrapolated  LiuVmann  process,'  ac¬ 
cording  to  which  the  residues  £.#  and  as  computed  from 

J  « 


Eqs.  (17)  and  (18)  are  regarded  as  a  measure  of  the 
"error"  of  the  estimated  solution,  and  are  "corrected 
by  an  amount  depending  upon  these  residues  according 
to  the  process. 


where  the  subscripts  denote  successive  approximations 
in  the  relaxation  process.  Here  is  a  coefficien 
of  over -relaxation,  defined  as, 

*=  £(«F  +  *lO 


where**,,  is  the  over-relaxation  coefficient  of 
.’''unkel, '  designed  to  give  optimum  convergence, 


*e  *  1  f  1+  cs  3.  ^1 

1  l  1  t  \  J  , 

with  I  ,  equal  to  the  number  of 

grid  intervals  along  the  %  and  ^  axes,  respec¬ 
tively.  in  Eq.  (21)  is  the  so-called  Richardson 

relaxation  coefficient,  a  constant,  and  equal  to  1/4. 
For  the  grid  employed  in  the  present  solutions,  the 
average  of  these  two  coefficients  was  found  to  0.34, 
and  was  employed  throughout  the  calculations.  Use 
of  this  coefficient  undoubtedly  speeded  the  conver¬ 
gence  of  the  relaxation  relative  to  the  Rich-rdson 
process,  at  least  during  the  earlier  stages;  the 
use  of  the  full  Frankel  coefficient  o*f  was  not 
attempted,  although  it  is  possible  that  this  may 
have  speeded  the  relaxation  convergence  3till  more. 
Charney  and  Phillips-5  have,  in  face,  found  the  use 


of  such  an  average  over-relaxation  coefficient  to  prod¬ 
uce  the  most  satisfactory  convergence  in  the  solution 
of  similar  equations. 

The  relaxation  cycle,  as  represented  by 
Eqs.  (19)  and  (20),  was  repeated  in  successive  rov;- 
by-rov;  scans  of  the  grid  until  the  residues  everywhere 
and  simultaneously  satisfied  the  condition, 

ISjjl  ^  €  » 

at  which  point  the  approximations  ^  / ft  and 

were  assumed  to  be  sufficiently  close  to  the 
actual  tendencies.  With  £  s  $4:  dUty 
representing  the  maximum  disposable  residue,  an 
average  of  six  scans  of  the  grid  -jere  required  to 
satisfy  this  convergence  condition  at  each  stage  of 
the  solution. 

Upon  completion  of  the  relaxation  as 
described  above,  the  iteration  process,  or  the  ad¬ 
vance  of  the  solutions  in  time,  was  carried  out  by 
use  cf  the  "centered"  extrapolation  formula. 


'lr'v'+l  =.  v'*-'  +  2  M  {ft)'*'  (  U3) 


..here  Avc  I,  j,  *  denotes  the  iterat  ion 

stage-,  the  increment  of  time  over  which  the 

solutions  are  to  be  extrapolated,  and 
the  variables  of  the  thermotropic  model .  Having 
selected  the  grid  point  spacing  -forehand 

to  faithfully  portray  the  large-scale  disturbances, 
the  time  increment  Ai  was  selected  one  hour 
in  order  to  satisfy  the  Courant-FrledKichss-Lo.’.y 
criterion  of  computational  stability.  For  two- 
parameter  or  two-layer  models  of  the  ; orl  hero 
considered,  and  for  the  iteration  scheme  of  Kq. t?>) * 
the  condition  that  the  solutions  not  gr.  ..  exponen- 
ent tally  with  time  is  approximately  of  the  form, 

where  C  is  to  be  interpreted  no.  the  maxi¬ 

mum  particle  speed  In  the  flow.  >,i3  Tb.e  selection 
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of  Uf  proved  computationally  stable  for  all 
cases  examined.  (See,  however,  footnote,  page  18.) 

For  the  first  hour’s  forecasts,  Eq.  (23) 
must,  of  course,  be  replaced  by  a  "forward"  extra¬ 
polation  formula. 


v*  =  nr 


6 


+  A  t 


(  \c 


(2.4) 


The  lateral  boundary  conditions  required 
in  the  solution  of  Eqs .  (9)  and  (10)  were  taken  as 
the  observed  24-hour  changes  centered  at  the  initial 
instant  ■fc-so  in  an  effort  to  specify  the  boundary 
tendencies  during  the  initial  stages  of  the  solution 
as  correctly  as  possible.  This  is  probably  an 
adequate  boundary  condition  for  some  research  pur¬ 
poses.  In  operational  application,  however,  a 
boundary  tendency  obtained  solely  from  past  data, 
a  synoptically -estimated  tendency,  or  a  zero  tend¬ 
ency  will  have  to  be  used  for  this  type  of  solution. 
To  permit  the  regeneration  of  initial  &  and  A 
data  on  the  original  18  x  23  grid,  these  boundary 
conditions  were  specified  on  the  three  outermost 
grid  rows,  and  were  held  fixed  throughout  the 
forecast  period  of  24  hours. 

2  3  The  Program  for  Machine  Solution 

The  finite-difference  equations  for  the 
thermotropic  model  were  programmed  and  coded  for 
the  IBM  Model  701  electronic  calculator,  a  high 
speed  digital  computer  of  the  stored  program  type.* 
The  program  itself  was  divided  into  several  log¬ 
ical  blocks  or  subroutines  as  follows: 

(1)  Basic  initialization.  This 
portion  of  the  program  prepares  and  positions  the 
initial  data  input,  and  sets  various  program 
controls  and  internal  tests  Into  readiness. 


*  The  basic  engineering  and  performance  charac¬ 
teristics  of  the  701  computer  are  described 
and  compared  with  those  of  similar  machines  in 
Proc.  J.R.E.,  Vol .  41,  No.  10,  Oct.  1953. 
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\<*))  Calculation  of  non-homogGuoouu  terms. 
This  .sab -program  calculates  the  several  members  of 
^  and  In  Kqu.  {<))  and  (10)  usints  the 

f liuto-dif Porenoou  dinnujiiud  errller . 

(3)  Calculation  of  x'eciduea  nnd  rela;;.- 
it'.on.  This  portion  of  the  program  computes  the 
rot,  1C  uni  a  from  F?qs: .  (1 '))  eat:  (1&)  corresponding  to 
the  variables  a  and  h,  us:’ ng  the  initially-loaded 
obiijfvvod  changes  ao  a  fii  ct  approximation  for  t  0 
us  a  tin.*  oolut  .onb  at  far  i= n*  At  •  The 

relaxation  is  then  performed  accordin’/  to  Eqs,  (19) 
•Mid  v20),  tho  solution  of  each  equation  being 
eomoietod  separately. 


(4)  Iteration, 
v^ut  tno  extrapolations  of 
taorc-hy  represent s from 
actual  numerical  "integr-* 
equation!.'  in  time . 


This  subprogram  carries 
Eqs.  'v2p )  arid  (24)  and 
.'iic  point  of  view,  the 
t'or.y'*  of  the  •overling 


(b)  Calculation  of  vertical  motion. 

Thiw  portion,  of  the  program  ;  as  utilized  only  at 
t  -  0,  12,  24  hours,  and  computes  the  average  vor¬ 
tical  velocity  according  to  Eq.  '16).  These  com¬ 
putation,':.  at  t  =  0  are  bhown  irA  the  Aopeudix . 


Each  of  these  subroutines  is  character¬ 
ized  by:  (a)  an  initial  series  of  orders  or 
machine  Instructions  in  which  the  addresses  (data 
locations  in  the  computer's  electrostatic  memory) 
of  the  pertinent  data  and  constants  are 
'' initialized"  or  set  to  the  beginning  of  the  finite - 
difference  grid;  (b)  a  sequence  of  instructions  by 
means  of  which  the  actual  calculations  is  performed — 
e.g.  addition,  multiplication  and  (c)  a  series  of 
instructions  v/hich  "modify"  or  alter  the  addresses 
of  the  data  used  in  the  calculations  in  order  that 
the  computation  may  proceed  to  the  next  grid  point.. 
The  sequence  of  steps  (b)  and  (c)  is  repeated 
cyclically  until  the  required  computation  has  been 
performed  for  each  point  of  the  finite -difference 
grid . 


While  the  approximate  character  of  the 
results  from  the  machine  integration  of  the  ther¬ 
motropic  model  could,  of  course,  be  anticipated 
to  some  extent,  It  was  found  necessary  to  carry  out 
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a  careful  manual  calculation  of  the  instantaneous 
tendencies  '«>%  /3t  and  in  order  to  obtain  a 

detailed  verification  of  the  correctness  of  the 
design  and  operation  of  the  code.  These  calcula¬ 
tions  allowed,  furthermore,  an  examination  of  the 
round-off  or  truncation  error  allowable  in  the 
machine,  and  made  possible  a  more  nearly  optimum 
selection  of  the  constants  oc  and  fi  and  a 
consequent  speeding  of  the  relaxation  process. 

While  executing  the  thermotropic  code 
for  a  24 -hour  forecast  at  500  and  lOOOmb,  the 
machine  performs  about  500,000  multiplications, 
approximately  half  of  which  are  made  in  the  course 
of  the  relaxation  process,  performed  each  hour  for 
Eqs .  ;’9)  and  (10).  In  addition,  the  machine  exe¬ 
cutes  approximately  15,000,000  additional  orders 
in  the  course  of  the  forecasts.  All  of  the  cal¬ 
culations  are  performed  completely  in  the  machine's 
high-speed  electrostatic  memory,  without,  the 
assistance  of  auxiliary  storage  units  such  as 
magnetic  drums  or  magnetic  tapes;  for  a  larger 
program  or  when  a  larger  finite-difference  grid  is 
involved,  the  use  of  these  lower-speed  storage 
units  will  probably  be  a  necessary  Ingredient  of 
ihe  program .  Ir,  its  present  form,  the  thermotropic 
program  requires  approximately  25  minutes  for  the 
production  of  12 -and  24-hour  forecasts  at  500  and 
1000  mb,  together  with  the  forecast  vertical  mo¬ 
tions;  this  time  includes  about  one  minute  for 
loading  the  program  into  the  calculator,  and 
approximately  three  minuter  for  the  automatic 
printing  of  the  forecasts.  In  all,  a  total  of 
sixty  12 -hour  forecasts  and  sixty  24 -hour  forecasts 
•were  prepared  at  both  500  and  1000  mb  from  this 
code,  starting  with  the  1500  G.C.T.  synoptic  chart 
for  i  Jan.,  1953,  and  including  each  0500  and  1500 
G.C.T.  chart  through  0500  G.C.T.  51  Jan.,  1953. 

2.*  Integration  of  the  Darotropic  Model 

As  the  second  model  in  this  series  of 
comparative  numerical  integrations,  the  barotropic 
model  representing  the  behavior  of  a  frictionless 
adiabatic  atmosphere  In  two-dimensional,  nor.- 
divergcnl  motion  over  flat  terrain  under  hydrostatic 
and  quasi -geostrcphlc  balance  was  examined.  This 


14 


model  may,  of  course,  be  recognized  as  a  special  case 
of  thermotropic  flow.  When  the  500  mb  surface  is 
identified  with  an  equivalence  level  representing  the 
mean  motion  of  the  atmosphere,  the  governing  equation 
of  the  barotropic  model  may  be  written. 


where  the  symbols  are  as  defined  earlier.  For  the 
barotropic  model,  the  equivalence  factor  c  was  taken 
as  unity.  Upon  the  introduction  of  a  scheme  of  non- 
dimensional  measure  and  finite-differences  as  for 
the  thermotropic  model,  Eq.  (25)  may  be  written. 


The  solution  of  Eq.  (2b)  proceeded  exactly 
as  in  the  case  of  the  thermotropic  equation  for  500 
mb  flow.  After  the  computation  of  residues,  the 
relaxation  proceeded  with  the  same  lateral  boundary 
conditions,  over-relaxation  coefficient  OC  , 
maximum  disposable  residue  6  ,  and  the  same 

iter?.*' ion  process  as  for  the  thermotropic  solutions. 
The  code  of  approximately  600  instructions  was 
written  by  Major  Zartner  and  requires  an  average  of 
five  minutes  for  execution  on  the  IBM  701  computer, 
including  the  performance  of  the  leading  and  print¬ 
ing  routines.  In  the  barotropic  case  there  is,  of 
course,  only  a  forecast  of  500  mb  height  changes, 
and  no  information  regarding  the  1000  mb  changes 
and  associated  vertical  motion  was  obtained. 

3.  The  Statistical  Summary  and  Analysis  of  Forecasts 

3.1  Introduction 

Upon  completion  of  the  integration  of  the 
barotropic  and  thermotrop!..  models  for  the  series 
of  60  cases,  it  became  a  primary  concern  to  es¬ 
tablish  some  measure  of  the  accuracy  of  the  thermo¬ 
tropic  and  barotropic  forecasts,  and  to  compare 
the  results  of  integrating  the  thermotropic  equations 
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with  those  obtained  from  the  barotropic  equation. 
Although  meteorologists  generally  recognize  the  short¬ 
comings  of  the  linear  correlation  as  a  measure  of 
'  forecasting  accuracy  it  is  nevertheless  the  most 

commonly  employed  statistic.  It  was  therefore  de¬ 
cided  to  compute  the  linear  correlation  between  the 
forecast  and  observed  height  changes  for  each  series 
of  numerical  forecasts.  Since  the  correlation  is 
principally  sensitive  to  phase  differences  the  root 
mean  square  error  was  computed  in  addition  to  fur¬ 
nish  a  more  amplitude -sensitive  measure  of  accuracy. 

It  is  unfortunate  that  the  correlations 
obtained  here  are  not  directly  comparable  with  those 
presented  by  other  investigators,  since  the  results 
depend  on  the  size  and  location  of  the  forecast  region 
and  on  the  nature  of  the  boundary  conditions  used  for 
a  particular  calculation.  In  fact,  the  possibility 
of  making  such  a  controlled  comparison  was  a  major 
motivation  of  the  present  work  in  which  the  solutions 
of  two  numerical  prediction  models,  the  thermotropic 
and  barotropic,  were  produced  under  identical  con¬ 
ditions  . 

While  Charney  and  Phillips^  have  com¬ 
pared  correlations  for  a  barotropic  model  with  those 
for  a  "2-  1/2  dimensional"  model,  the  barotropic 
model  produced  a  forecast  at  the  500  mb  level  and 
the  2-  1/2  dimensional  model  produced  forecasts  at 
i  the  pOO  and  7’Q0  mb  levels.  The  quantities  compared 

j  were  the  correlations  of  the  barotropic  forecast 

500  mb  changes  with  the  observed  500  mb  changes  and 
the  correlations  of  the  2-  1/2  dimensional  model 
forecast  700  mb  changes  with  the  observed  700  mb 
changes.  Charney^  has  more  recently  presented  a 
comparison  of  the  500  mb  forecasts  of  a  three-level 
model  with  500  mb  forecasts  of  a  one-level  model  and 
the  700  mb  forecasts  of  a  two-level  model  during  a 
three-day  period  of  Intense  cyclogenesis.  The 
correlations  of  forecast  24-hour  changes  with  ob¬ 
served  24 -hour  changes  were  high  for  the  three  series 
of  forecasts  preceding  the  storm,  but  only  those 
for  the  three-level  forecasts  remained  high  after 
the  storm  developed.  Neither  of  these  compar¬ 
isons  were  felt  to  constitute  a  conclusive 
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examination  of  the  relative  performance  of  the  one-, 
two-,  and  three-level  models,  although  they  docu¬ 
mented  for  the  first  time  the  anticipated  superior¬ 
ity  of  the  three-level  model  in  a  case  of  strong 
baroclinity.  There  remains  to  be  made,  however, 
systematic  model  comparisons  in  a  wide  variety  of 
synoptic  situations. 

From  the  present  work  we  are  in  a  position 
to  compare  the  performances  of  the  barotropic  and 
thermotropic  models  in  predicting  500  mb  height 
changes  over  a  long  series  of  forecasts,  in  which 
both  models  were  integrated  under  identical  con¬ 
ditions  . 


In  the  following  sections,  the  several 
basic  statistics  (correlations,  root -mean -square 
errors)  summarizing  the  results  of  the  integrations 
will  be  presented  and  briefly  described.  Further 
statistical  tests  conducted  will  be  discussed; 
finally,  an  evaluation  of  the  "skill"  of  the  fore¬ 
casts  from  a  statistical  viewpoint  will  be  made. 

5.2  Space  Statistics 

Linear  correlations  between  the  forecast 
and  observed  height  changes  and  the  root  mean 
square  errors  of  the  forecast  height  changes  were 
computed  for  the  12-  and  24 -hour  500  mb  barotropic 
and  thermotropic  forecasts,  and  the  12-  and  24- 
hour  1000  mb  thermotropic  forecasts.  In  addition 
to  the  space  correlation  and  space  root  mean  square 
error  for  each  forecast,  a  time  correlation  and 
time  root  mean  square  error  at  each  grid  point  were 
computed  for  each  series  of  forecasts.  All  of  the 
computations  were  performed  on  the  IBM  type  701 
electronic  computer. 

Figures  3.1  through  3.6  are  the  space 
correlations  for  each  series  of  forecasts  plotted 
as  a  function  of  time;  Figures  3.7  through  3.12 
are  the  space  root  mean  square  errors  plotted  as 
a  function  of  time.  The  most  striking  feature  of 
the  time  distribution  of  the  space  correlations  is 
to  be  noted  by  comparing  the  results  for  the  baro¬ 
tropic  and  thermotropic  500  mb  predictions  (Figures 
3.1  through  3.4).  For  the  12-  and  24 -hour 
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predictions  the  barotropic  and  thermotropic  space  cor¬ 
relations  are  almost  identical.  In  those  cases  for 
which  the  space  correlation  of  the  barotropic  pre¬ 
diction  is  low  the  space  correlation  of  the  thermo¬ 
tropic  prediction  is  also  low.*  Two  possible  ex¬ 
planations  may  be  offered  for  this  result: 

(1)  Errors  in  the  prediction  may  be  due  to 
the  failure  to  include  some  pertinent  physical  process, 
in  which  case  it  must  be  concluded  that  the  omission 
occurs  in  both  models. 

(2)  There  is  some  overpowering  influence, 
which  may  not  be  a  physical  influence,  that  causes 
both  models  to  behave  similarly.  This  may  be  an 
effect  due  to  the  mathematical  formulation  of  the 
problem,  or  due  to  the  use  of  similar  numerical 
approximations  in  the  solutions  of  the  tw^  models. 


One  exception  to  this  behavior  is  to  be  noted  in 
the  24 -hour  500  mb  prediction  made  at  OJOOZ  17 
January.  The  space  correlation  is  much  lower  for 
the  bar  ’Topic  forecast  than  for  the  thermotropic 
forecast;  there  is  a  corresponding  difference  in 
the  space  root  mean  square  errors.  This  difference 
is  due  to  a  very  large  error  in  the  24 -hour  baro¬ 
tropic  forecast,  which  is  concentrated  in  the 
northeast  corner  of  the  forecast  area.  This  error 
is  as  yet  unexplained,  but  since  it  arose  in  the 
computations  after  a  successful  12-hcur  forecast 
had  been  made,  it  is  conceivable  that  the  error 
is  a  manifestation  of  some  form  of  computational 
Instability . 
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The  general  level  of  the  24 -hour  500  mb  space 
correlations  (Figures  5-3  and  3.4)  is  only  slightly 
lower  than  the  12-hour  500  mb  space  correlations 
.fFigures  5.1  and  3.2)  except  for  the  forecasts  made  at 
15Z-7  Jan.,  15Z-18  Jan.,  03Z-19  Jan.,  15Z-19  Jan.,  and 
15Z-28  Jan.,  when  the  24-hour  correlation  is  much 
lower  than  the  12 -hour  correlation.  This  implies  that 
the  forecast  calculations  compound  an  error  in  the 
change  pattern.  A  forecast  which  starts  off  badly 
during  the  first  12  hours  goes  further  astray  during 
the  next  12  hours;  a  forecast  which  starts  off  well 
during  the  first  12-hour  period  will  be  a  good  fore¬ 
cast  at  the  end  of  the  24-hour  period.  A  comparison 
of  the  graphs  of  the  space  RMSE  for  the  barotropic 
and  thermotropic  500  mb  predictions  (Figures  3.7 
through  3.10)  reveals  the  same  striking  similarity 
present  in  the  graphs  of  the  correlation  coefficients. 
The  RMSE  for  the  24 -hour  forecasts  are  greater  than 
those  for  the  12-hour  forecasts  in  all  cases.  The 
correspondence  between  cases  for  which  the  RMSE  of 
both  the  12-  and  24-hour  predictions  are  relatively 
large  may  be  noted . 

The  12 -hour  500  mb  barotropic  and  thermo¬ 
tropic  space  correlation  coefficients  were  examined 
in  some  detail.  The  average  correlation  coefficient 
for  the  barotropic  predictions,  0.817,  was  slightly 
higher  than  the  value  for  the  thermotropic  predic¬ 
tions,  0.312;  the  difference  between  these  correla¬ 
tions  is  not  statistically  significant.  The  differ¬ 
ences  between  the  daily  values  of  the  normally  dis¬ 
tributed  Fisher's  Z  for  the  barotropic  and  for  the 
thermotropic  predictions  were  compared  to  an  esti¬ 
mate  of  the  standard  error  based  on  a  population  of 
204  (number  of  grid  points)  independent  pairs. 

Since  the  number  of  Independent  pairs  is  considerably 
less  than  204  this  procedure  leads  to  an  underestimate 
of  the  standard  error.  The  daily  differences  of  the 
Z  values  exceeded  twice  the  estimated  standard  error 
in  only  three  of  the  60  cases  (5#) .  Since  the 
standard  error  was  considerably  underestimated  it 
must  be  concluded  that  the  barotropic  and  thermo¬ 
tropic  3pace  correlation  coefficients  for  the 
12-hour  500  mb  forecasts  are  not  significantly 
different  for  any  of  the  60  predictions. 
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The  space  correlation  coefficients  for  the 
1000  mb  thermotropic  predictions  are  generally  lower 
than  those  for  the  500  mb  predictions.  The  general 
level  of  the  1000  mb  space  correlations  is  about  the 
same  for  the  12-  as  for  the  24-hour  predictions, 
except  that  the  minimum  24 -hour  correlations  are 
considerably  lower  than  the  minimum  12 -hour  corre¬ 
lations  (Figures  3-5  and  3-6).  The  correspondence 
between  poor  12 -hour  predictions  and  poor  24 -hour 
predictions  is  not  as  well  marked  as  for  the  500  mb 
forecasts  .  The  space  root  mean  square  e  '•s  of  the 
24-hour  1000  mb  predictions  are  greater  than  those 
for  the  12-hour  predictions  for  all  cases;  the  in¬ 
crease  in  error  from  the  12-  to  24-hour  forecast  is 
almost  linear  in  time. 

The  frequency  distributions  of  the  space 
correlation  coefficients  are  presented  in  Figures 
3.13  through  3.18.  A  comparison  of  the  frequency 
distributions  of  the  space  correlations  for  the  500 
mb  thermotropic  and  barotropic  forecasts  (Figures 
3.13  through  3.16)  again  illustrates  the  similarity 
of  the  two  sets  of  forecasts.  The  distributions  of 
the  1000  mb  space  correlations  are  flatter  than  the 
distributions  of  the  500  mb  space  correlations. 

The  median  value  of  the  space  correlation  for  each 
set  of  forecasts  is  tabulated  in  Table  3.1.  From 
each  frequency  distribution  the  percentage  of  fore¬ 
casts  for  which  the  space  correlation  exceeded  0.7 
may  be  computed.  In  these  cases  the  forecast 
method  explained  approximately  5*#  or  more  of  the 
variance  of  the  change  pattern.  These  percentages 
are  also  given  in  Table  1.  Considerably  mor<.  th  n 
half  of  the  barotropic  and  thermotropic  50Q  mb 
forecasts  exceed  the  0.7  correlation  level:  loss 
than  half  of  the  thermotropic  1000  mb  forecasts 
exceeded  the  0.7  correlation  level. 

In  concluding  our  discussion  of  the 
space  correlations  3ome  general  observations  may  be 
noted.  The  result  that  the  correlation  is  higher 
if  the  actual  height  changes  are  larger  (in  absolute 
magnitude)  has  been  observed  and  reported  by 
Lonnqvlst  10  in  the  case  of  standard  forecas*-  tech¬ 
niques  ami  by  the  Staff  Members,  Institute  of  Meteoro¬ 
logy,  University  of  Stockholm11  in  the  case  of 
numerical  forecasts  using  the  barotropic  model. 
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Table  3  .1  --  Space  correlations.  Median  values  and 
percentage  of  cases  for  which  corre¬ 
lation  was  greater  than  0.7. 


Level Model  Period  Median  Value  Percentage 
_ Exceeding  0.7 

12  Hour  .82  75.0 

Thermo¬ 

tropic 

24  Hour  .80  75-3 

500  mb 

12  Hour  .83  78.3 

Baro- 

tropic 

24  Hour  .80  66.6 


12  Hour  .68  46.7 


Thermo - 

1000  mb  tropic 


24  Hour 


.69 


48.3 


Table  3.2  --  Average  space  root -mean-square  errors. 


Level 


500  mb 


1000  mb 


Model 

Period 

Mean  RMSE 
in  Feet 

Thermo¬ 

12  Hour 

118.7 

tropic 

24  Hour  . 

231.2 

Baro- 

12  Hour 

115.8 

troplc 

24  Hour 

228.1 

Thermo¬ 

12  Hour 

126.4 

tropic 

24  Hour 

223.1 
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The  same  general  trend  has  been  observed  in  the  pres¬ 
ent  investigation.  This  is  probably  a  manifestation 
of  the  tendency  to  obtain  a  high  correlation  between 
two  populations  if  each  population  has  a  bimodal  dis¬ 
tribution.  If  the  average  absolute  value  of  the 
height  change  for  a  particular  case  is  high,  the 
field  is  generally  one  of  relatively  large  rises  and 
falls  of  contour  height,  with  only  a  small  area  of 
little  change.  The  nature  of  the  linear  correlation 
is  such  that  a  higher  correlation  is  to  be  expected 
in  such  cases  than  in  cases  for  which  the  height 
change  field  is  of  smaller  contrast. 

The  average  space  root  mean  square  errors 
are  given  in  Table  3.2.  We  again  note  the  similarity 
between  the  results  of  the  500  mb  thermotropic  and 
barotropic  predictions.  The  average  space  root  mean 
square  errors  for  the  12-  and  24 -hour  1000  mb  fore¬ 
casts  are  approximately  equal  to  the  averages  for 
the  corresponding  500  mb  forecasts.  This  means  that 
the  1000  mb  predictions  are  somewhat  poorer  than  the 
500  mb  predictions  since  the  standard  deviation  of 
the  observed  changes  of  the  height  of  the  500  mb 
surface  is  approximately  30#  larger,  in  the  average, 
than  the  standard  deviation  of  the  1000  mb  changes. 

This  is  further  Indicated  by  the  result  that  the 
space  correlations  for  the  1000  mb  predicted  changes 
are  considerably  lower  than  those  for  the  5°0  mb 
predicted  changes.  The  average  RMSE  for  the  24- 
hour  forecasts  is  approximately  twice  the  average 
RMSE  for  the  corresponding  12-hour  forecasts. 

The  error  in  magnitude  of  the  forecast  changes  is 
essentially  linear  in  time. 

3.3  Time  Statistics 

The  geographical  distributions  of  the  time 
correlation  coefficients  are  presented  in  Figures 
3.19  through  3.24.  For  each  series  of  forecasts  the 
forecast  changes  at  a  grid  point,  ordered  sequentially 
in  time,  were  correlated  with  the  observed  changes, 
in  the  same  order,  at  the  same  grid  point;  this  yields 
correlation  coefficients  for  all  of  the  grid  points 
as  a  measure  of  the  geographical  variation  of  the 
forecast  accuracy.  There  are  some  features  of  the 
geographical  pattern  common  to  each  forecast  series. 

Each  figure  has  two  pronounced  maxima,  one  over  the 
eastern  U.  S.,  and  the  other  over  the  southwestern  U.  S. 
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Fig.  3.24 
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There  is  a  minimum  value  over  the  northwestern  part 
of  the  forecast  area  which  entends  southeastward 
into  the  central  U.  S.  The  distributions  of  the 
time  correlations  for  the  barotropic  and  thermo¬ 
tropic  12 -hour  500  mb  predictions  (Figures  3.19 
and  3.20)  are  strikingly  similar  (as  were  the  time 
distributions  of  the  space  correlations) .  A 
comparison  of  these  sets  of  time  correlations  indi 
cates  that  there  is  no  grid  point  for  which  the 
difference  between  the  correlations  is  statisti¬ 
cally  significant.  The  chart  of  the  time  correla¬ 
tions  for  the  24-hour  500  mb  barotropic  predic¬ 
tions  (Figure  3.21)  differs  from  the  chart  of  the 
time  correlations  for  the  24-hour  500  mb  thermo¬ 
tropic  predictions  in  the  northeast  corner  of  the 
forecast  area  where  the  barotropic  correlations 
have  a  minimum.  This  minimum  is  caused  by  one 
very  poor  barotropic  24-hour  prediction  (see  foot¬ 
note,  page 


A  comparison  of  the  charts  for  the  12- 
hour  500  mb  forecasts  and  the  24 -hour  500  mb 
forecasts  indicates  that  the  area  of  high  correla¬ 
tion  in  the  eastern  part  of  the  forecast  region  is 
smaller  for  the  24 -hour  forecasts  than  for  the  12- 
hour  forecasts.  Consider  the  area  with  correlation 
greater  than  0.85;  its  western  edge  is  at  approx¬ 
imately  the  same  position  for  the  12-  and  24 -hour 
300  mb  predictions,  but  Its  northern  extent  is 
diminished.  The  v;estern  area  of  correlation 
greater  than  O.85  is  considerably  smaller  for  the 
24 -hour  500  mb  forecasts  than  for  the  12-hour  500  mb 
forecasts.  The  minimum  value  of  the  correlation  is 
lower  for  the  24 -hour  predictions  than  for  the  12- 
hour  predictions.  The  area  of  correlation  lower 
than  0.7  is  larger  on  the  charts  of  the  24-hour 
500  mb  time  correlations  than  on  the  charts  of  the 
12-hour  500  mb  time  correlations.  The  decrease 
of  the  correlation  from  12-  to  24-hours  is  espe¬ 
cially  marked  in  the  northwest  portion  of  the  fore¬ 
cast  area  and  in  the  "trough"  east  of  the  Kockies. 
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The  principal  features  of  the  pattern  of  the 
time  correlation  coefficients  for  the  1000  mb  thermo¬ 
tropic  forecasts  (Figures  3.23  and  3.24)  are  quite 
similar  to  the  patterns  of  the  500  mb  correlations. 

The  maximum  values  of  the  correlation  are  of  the  same 
magnitude  as  those  for  the  500  mb  predictions,  but 
the  minimum  values  are  much  lower.  The  result  is  that 
the  gradient  of  correlation  for  the  1000  mb  forecasts 
is  mucn  greater  than  on  the  charts  for  the  500  mb 
forecasts.  The  region  of  low  correlation  in  the  north¬ 
west,  extending  southeastward  into  the  region  east  of 
the  Rockies  in  the  U.  S.,  is  much  more  pronounced. 

The  general  result  that  the  correlation  of 
observed  changes  with  forecast  changes  is  higher  for 
the  eastern  part  of  the  forecast  region  than  for  the 
western  part  has  also  been  obtained  by  Bushby  and 
Hinds1  for  integrations  of  the  Sawyer-Bushby  two- 
parameter  baroclinic  model.  The  forecast  region  in 
thiis  case  covered  the  North  Atlantic  Ocean  and 
Western  Europe.  The  result  was  Indicated  for  both 
tne  500  and  1000  mb  forecasts,  although  only  a  limited 
number  of  forecasts  were  made.  Bushby  and  Hinds  showed 
that  the  space  correlations  were  higher  if  only  the 
eastern  portion  of  the  forecast  region  was  used  in 
the  verification  than  if  the  entire  forecast  region 
was  used.  The  Staff  Meqib  !s,  Institute  of  Meteorology, 
University  of  Stockholm1  ,.ave  presented  a  chart  of 
the  time  correlations  bt cween  observed  and  computed 
500  mb  height  changes  for  a  somewhat  longer  series 
of  barotropic  forecasts  over  approximately  the  same 
area,  and  find  the  same  increase  of  correlation 
coefficient  in  the  eastern  portion  of  the  forecast 
region.  In  both  of  these  series  of  forecasts 
'operational"  assumpt ions  were  made  concerning  the 
boundary  conditions  so  that  only  in  ormation  avail¬ 
able  at  the  start  of  the  forecast  period  was  used. 

The  geographical  distributions  of  the  time 
root -mean-3quare  errors  are  shown  in  Figures  3.25 
through  3.30.  The  patterns  are  partiallv  obscured 
by  a  north-south  gradient  introduced  by  "climatology." 
Since  the  variability  of  the  height  of  a  constant 
pressure  surface  increases  toward  the  north  of  the 
forecast  area,  the  magnitude  of  the  error  in  the 
forecast  height  change  Increases  toward  the  north. 

This  effect  superimposes  a  zonal  pattern  on  the  root- 
mean-sq  ;are  error  charts.  It  is  still  possible  to 
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distinguish  some  features  apart  from  the  zonal  pattern. 
The  errors  along  the  northeast  coast  of  the  U.  S.  are 
somewhat  smaller  than  those  along  the  Pacific  north¬ 
west  for  the  12-hour  500  mb  predictions;  the  maximum 
errors  in  eastern  Canada  are  smaller  than  the  maximum 
errors  in  western  Canada.  A  comparison  of  the  charts 
for  the  12-  and  24-hour  500  mb  predictions  shows  the 
same  linearity  in  time  of  the  root-mean-square  error 
that  was  noted  in  the  space  root -mean-square  errors. 

The  space  distribution  of  the  time  root- 
mean-square  errors  of  the  24 -hour  500  mb  thermotropic 
predictions  (Figure  3.28)  shows  that  the  error  along 
most  of  the  eastern  edge  of  the  forecast  area  is  less 
than  the  e^ror  along  the  western  edge.  The  maximum 
error  over  _astern  Canada  is  less  than  the  maximum 
error  over  western  Canada.  The  geographical  distri¬ 
bution  of  the  time  root -mean-square  errors  of  the 
24-hour  500  mb  barotropic  predictions  shows  a  pro¬ 
nounced  maximum  at  uhe  northeast  corner  of  the  fore¬ 
cast  region.  These  large  errors  are  caused  by  one 
very  poor  barotropic  forecast  (see  footnote,  page 

The  geographical  distributions  of  the  time 
root -mean-square  errors  of  the  1000  mb  thermotropic 
forecasts  correspond  strongly  to  the  time  correlation 
patterns,  expeclally  in  the  areas  of  poor  forecast 
(cf.  Figures  3.23  and  3.29,  3.24  and  3.30).  The 
locations  of  the  time  root -mean-square  error  maxima 
are  in  the  same  region  as  the  minima  of  the  time 
correlations.  This  region  is  the  principal  feature 
of  the  1000  mb  time  root -mean-square  error  patterns. 
Each  chart  also  shows  regions  of  error  minima  at  the 
southern  corners  of  the  forecast  area. 

3.4  Normalized  Time  Root -Mean-Square  Error 

The  time  root -mean-square  errors  may  be 
normalized  in  order  to  remove  the  bias  toward  a  purely 
zonal  pattern.  It  shall  first  be  hypothesized  that 
the  standard  deviation  of  the  time  RMSE  in  one  row 
of  grid  points  is  approximately  equal  to  the  standard 
deviation  that  would  be  obtained  for  each  grid  point 
in  the  row  if  values  of  the  RMSE  were  available  from 
many  sets  of  60  forecasts.  It  is  also  assumed  that 
the  values  of  the  RMSE  are  normally  distributed. 

The  normalized  RMSE  shall  then  be  defined  by: 
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where  RMSE  is  the  mean  of  ail  values  in  a  row  and 
O'  (RMSE)  is  the  standard  deviation  of  the  values 
in  the  row.  The  normalized  RMSE  is  such  that 
values  at  different  latitudes  may  now  be  directly 
compared  since  the  bias  in  the  RMSE  due  to  the 
increase  in  magnitude  of  the  observed  height  changes 
with  increasing  latitude  has  been  removed  .  A  point 
at  which  the  value  of  FMSS..  :i c  greater  than  one 
• RMSEfl  >  1)  is  a  point  for  .mien  the  forecasts  , ere 
poor  in  tne  sense  that,  the  value  of  the  FMSE  is 
greater  than  the  mean  FMSE  by  more  than  one 
standard  deviation.  A  difference  of  this  magnitude 
would  occur  randomly  in  less  than  33  cases  in  a 
hundred.  (Under  our  assumptions  this  means.  S3 
values  out  of  one  hundred  values  of  the  RMSE  com¬ 
puted  for  eacli  of  100  sets  of  60  forecasts  each.) 

A  value  of  the  RMSE.,  greater  than  two  is  a  "poor" 
value  in  the  sense  that  a  difference  of  this  magni¬ 
tude  would  occur  randomly  in  less  than  ; ! vc  cases 
in  a  hundred.  Values  of  the  P.MSE^  less  than  -1 
and  -2  arc  "good,"  In  the  same  sense,  at  the  ?■)% 
and  level,  respectively.  The  normalized  time 
RMSE  were  computed  for  the  <?4 -hour  thermotropic 
300  and  1000  mb  forecasts.  Values  of  the  normalized 
RMSE  In  different  geographical  regions  are  directly 
comparable  and  it.  Is  possible  to  examine  the  geo¬ 
graphical  pattern  of  the  normalized  FMSE  in  some 
tie  tail.  For  the  purpose  of  comparing  this  pattern 
with  the  eanh's  topography  in  the  forecast  re  .'ion 
a  special  base  map  with  contour  lines  drawn  to 
fit  values  of  the  terrain  elevation  averaged  over 
five  degree  latitude-longtitude  areas  was  prepared . 

The  patterns  cf  normalized  RMSE  for  the 
2 -i  -hour  thermotropic  1000  and  500  mb  forecasts  are 
presented  superimposed  on  the  topography  char*  in 
Figures  3*3 1  and  3.32,  respectively.  The  must 
striking  feature  is  the  location  of  the  region  of 
KMSEkj>  i  on  the  chart  for  1000  mb  (Figure  3.31)  in 
the  lee  of  the  Rocky  Mountains.  In  terms  cf  the 
measure  provided  by  the  normalized  RMSE  the  region 
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of  poorest  forecast  is  in  western  Kansas  south  to 
the  Texas  panhandle.  In  the  chart  for  500  mb 
(Figure  3.32)  the  region  of  RMSEN  >1  covers  the 
northern  half  of  the  western  boundary,  probably 
due  to  errors  introduced  at  the  boundary,  and 
continues  southeastward  in  the  lee  of  the  Rocky 
Mountains.  In  both  charts  the  forecasts  are  better 
in  the  eastern  half  of  the  region  than  in  the 
western  half,  except  for  the  area  of  RMSE,.  < -1  in 
the  southwestern  corner.  The  chart  for  trie  500  mb 
RMSEm  has  one  small  area  of  very  good  forecasts, 
(RMSEn  <  -2)  in  the  center  of  the  eastern  boundary. 
This  region  may  be  due  to  the  favorable  influence 
of  the  boundary  conditions. 

The  values  of  the  RMSE  in  a  particular 
column  of  grid  points  may  be  averaged  to  obtain 
a  mean  measure  of  the  west-to-east  variation  of 
th<=>  forecast  accuracy.  The  average  values  of 
RMSEjj  for  each  column  are  presented  in  Figure  3.33> 
where  the  mean  height  of  the  terrain  is  indicated 
by  the  shaded  area.  The  mean  RMSE  for  both  the 
500  and  1000  mb  forecasts  show  maxima  in  the  lee 
of  the  Rockies.  The  500  mb  curve  also  has  a 
secondary  maximum  in  the  lee  of  the  lower  Appalachian 
Range.  The  curves  presented  here  support  the  deduc¬ 
tion,  made  from  the  geogx’aphlcal  distribution  of 
the  time  correlations,  that  the  forecasts  are 
poorer  In  the  western  half  of  the  region  than  in 
the  eastern  half.  The  geographical  distribution 
of  the  normalized  time  root -mean-square  errors 
Indicate,  much  more  clearly  than  any  of  the 
other  data,  that  the  region  of  poorest  forecast  is 
in  the  lee  of  the  Rocky  Mountains.  This  result 
may  be  due  to  some  deficiency  in  the  numerical 
prediction  model.  There  are  two  classes  of 
physical  processes  omitted  from  the  prediction 
models  tthat  may  be  responsible.  One  may  be  termed 
a  purely  lee  effect;  the  surface  of  the  earth  has 
been  assumed  perfectly  flat  in  the  derivation  of 
the  prediction  equations.  No  consideration  has 
been  taken  of  the  fact  that  air  is  forced  to  flow 
up  over  the  mountains.  The  other  effect  is  one 
due  to  the  Rocky  Mountains  serving  as  a  crude 
physical  barrier,  channeling  the  low-level  outbreaks 
of  polar  air  along  their  eastern  slopes.  The 
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process  of  low-level  advection  is  not  included  in  the  pre¬ 
diction  model.  It  is  therefore  implied,  by  the  analysis  of 
the  present  results,  that  further  generalizations  of  the 
numerical  prediction  models  might  well  be  directed  toward 
the  inclusion  of  topographical  effects  and  low-level  ad¬ 
vection  processes. 

3.5  Comparison  With  No-Skill  Forecasts 

The  results  presented  up  to  this  point  leave  one 
question  unanswered:  "How  do  the  correlation  coefficients 
for  the  numerical  pretlicticns  compare  with  correlation  co¬ 
efficients  of  forecasts  obtained  by  some  other  technique?" 

In  seeking  such  a  comparison,  it  is  desirable  to  use  a 
"no-skill"  forecast  technique  which  utilizes  the  same  given 
information  as  the  numerical  forecasting  techniques.  The 
significance  of  the  comparison  may  be  increased  by  the  lack 
of  sophistication  of  the  no-skill  forecast  technique  and  the 
randomness  of  its  selection.  It  is  quite  probable  that,  if 
enough  different  schemes  are  tried,  we  might  eventually  dis¬ 
cover  a  simple  scheme  which  would  forecast  "as  well  as" 
numerical  prediction  (say,  in  terms  of  the  correlation  co¬ 
efficient  )  . 


In  view  of  these  considerations  .  e  shall  now  pro¬ 
pose  a  simple  forecast  technique.  The  first  step  in  preparing 
a  12-hour  forecast  is  to  shift  the  field  of  r'0Q  mb  height 
changes  for  the  12-hour  period  ending  at  the  start  of  the 
forecast  period  two  grid  intervals  to  the  east  (an  eastward 
velocity  of  approximately  7.5*  longitude  per  12  hours  at 
North) .  In  the  next  step  we  make  use  of  the  heights  at 
12-hour  verification  time  on  the  boundary  (in  order  to  make 
the  technique  comparable  with  the  numerical  forecasting 
technique  described  in  Section  2)  to  compute  the-  observed 
12 -hour  500  mb  boundary  height  changes.  The  boundary  changes 
u:v  tiun  used  to  Improve  he  interior  values  obtained  by 
straight  extrapolation.  An  Improved  forecast  12-hour  change 
for  the  first  interior  row  is  obtained  by  averaging  the  ex¬ 
trapolated  height,  change,  weighted  one,  and  the  observed 
height  change  at  the  neighboring  boundary  point,  weighted 
three.  The  value  for  the  second  row  is  obtained  by  averaging 
the  extrapolated  and  boundary  changes,  equally  weighted.  The 
third  interior  row  of  forecast  changes  is  obtained  by  aver¬ 
aging  the  extrapolated  change,  weighted  three,  and  the  bound¬ 
ary  change,  weighted  one.  In  each  corner  of  the  forecast 
region  there  are  nine  points  for  which  the  forecast  changes 
are  obtained  by  averaging  the  two  values  computed  from  each 
boundary.  The  scheme  outlined  provides  a  12 -hour  500  mb 
"forecast"  in  the  3ame  sense  that  the  numerical  predictions 
do.  The  interior  data  is  based  on  information  available  up 
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to  forecast  time;  boundary  data  is  available  at  verifi¬ 
cation  time . 

For  a  24-hour  forecast  no  further  boundary  data 
are  available;  the  forecast  height  change  pattern  for  the 
first  12-hour  period  is  shifted  eastward  two  grid  intervals 
as  a  forecast  for  the  second  12-hour  period.  The  forecasts 
for  the  first  and  second  12-hour  periods  are  added  to  yield 
the  24-hour  500  mb  height  change  prediction. 

The  "no-skill"  12-  and  24 -hour  500  mb  forecasts 
were  made  for  a  sample  of  30  randomly  selected  cases  from 
among  the  60  for  which  the  numerical  forecasts  were  made. 

T.ie  forecast  12-  and  24-hour  500  mb  height  changes  were  cor¬ 
related  with  the  observed  12-  ana  24-hour  500  mb  height 
changes,  respectively,  for  each  forecast.  The  root-mean- 
square  error  of  the  forecast  height  changes  was  also  computed. 
Average  values  of  these  statistics  are  compared  with  the 
averages  for  the  thermotropic  forecasts  for  the  same  30  cases 
in  Table  3-3.  On  a  daily  basis  the  thermotropic  2'+ -hour  500  mb 
forecasts  were  better  than  the  no-skill  forecasts  at  the  one 
percent  level  in  14  cases,  while  the  no-skill  forecasts  were 
better,  at  the  same  level,  in  only  three  cases.  The  frequency 
distributions  of  the  24-hour  space  correlations  for  the  thirty 
24 -hour  no-ski II  forecasts  and  the  corresponding  thermotropic 
forecasts  are  shown  in  Fig.  3*34.  The  no-skill  forecasts  did 
not  do  as  well  as  the  best  thermotropic  forecasts  nor  as 
poorly  as  the  worst,  i.e.  in  this  sense  it  is  a  more 


Table  3.3  --  Comparison  of  the  12-  and  24-hour  500  mb 
forecasts  obtained  by  the  no-skill  tech¬ 
nique  with  those  obtained  by  the  numerical 
integration  of  the  thermotropic  equations. 
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conservative  forecast  scheme. 

The  time  correlations  of  the  forecast  24 -hoar  500 
mb  height  changes  with  the  observed  height  changes  at  each 
grid  point  were  also  computed  and  compared  with  the  time 
correlations  of  the  24-hour  500  mb  thermotropic  forecasts. 
Figure  3.35  indicates  the  result  of  this  comparison.  In  the 
shaded  area,  located  to  the  east  of  the  center  of  the  fore¬ 
cast  region,  the  thermotropic  correlations  are  greater  than 
the  no-skill  correlations  by  an  amount  that  is  significant 
at  the  one  percent  level.  Near  the  boundary  there  is  no 
pronounced  superiority  of  one  scheme  over  the  other.  In 
the  interior  of  the  region,  where  the  no-skill  technique 
does  not  take  account  of  boundary  effects,  the  integration 
of  the  equations  is  a  much  better  technique  than  the  straight 
extrapolation  of  the  height  change  field.  There  are  two  pos¬ 
sible  explanations  for  this  result.  Either  the  equations  "do 
something"  in  the  Interior,  quite  independently  of  boundary 
influences, that  is  more  reliable  meteorologically  than  an 
extrapolation  of  the  past  change  pattern,  or  the  equations 
furnish  a  mechanism  for  propagating  the  boundary  information 
into  the  very  center  of  the  forecast  region. 

It  is  quite  possible  that  the  no-skill  forecasts 
might  be  improved  by  using  an  interpolation  scheme  that 
brings  • he  influence  of  the  observed  12-hour  changes  further 
iv. o  the  region.  Results  might  also  be  improved  if  observed 
f'00  mb  heights  at  the  12-hour  verification  time  were  used 
along  • hree  boundary  rows  (as  in  the  numerical  forecasts) 
instead  of  along  one  (as  in  the  no-skiil  forecasts).  In  this 
case  a  more  elegant  interpolation  scheme,  based  on  the 
gradient  of  the  observed  12-hour  change  normal  to  the  bound¬ 
ary.  might  be  used  to  improve  the  Interior  values. 

3.0  Summary 

The  correlations  of  the  forecast  500  mb  changes 
wi;h  observed  changes  are  high  enough  to  be  quite  encouraging; 
the  correlation  of  the  1000  mb  forecast  changes  with  observed 
changes  are  not  as  high  but  are  nevertheless  felt  to  be  sig¬ 
nificant  In  view  of  the  relatively  simple  physical  theory 
employed.  More  research  is  required  to  establish  clearly  the 
influence  cf  the  boundary  conditions  on  the  final  forecast. 

The  similarity  of  the  behavior  of  the  correlations  of  the 
SCO  mb  thermotropic  forecasts  and  the  correlations  of  the 
h 00  mb  barotropie  forecasts  remains  the  most  striking  feature 
of  the  results  presented  here.  It  is  to  be  hoped  that  further 
more  subtle,  analysis  of  the  numerical  forecasts  may  reveal 
the  distinctions  between  the  physical  behavior  of  the  two 
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models  which, have  been  indicated,  for  example,  bTr  stabil¬ 
ity  studies.  The  results  of  Section  3.4  indica  chat 
further  improvement  in  numerical  forecasting  may  be  gained 
within  the  framework  of  the  present  theory  by  including 
topographical  effects. 

4 .  The  Synoptic  Summary  and  Analysis  of  Forecasts 

4.1  Introduction 

To  date,  the  interpretation  of  most  numerical 
forecasts  has  been  from  the  point  of  view  of  the  dynamic 
meteorologist,  i.e.,  the  predictions  are,  in  fact,  solu¬ 
tions  to  a  set  of  hydrodynamical  equations  and  offer  direct 
evidence  of  the  accuracy  of  the  theory  and  the  suitability 
of  the  numerical  procedures  employed. 

In  addition  to  this  interpretation,  it  was  felt 
that  in  the  present  series  of  forecasts,  a  valuable  summary 
of  the  results  would  be  from  the  point  of  view  of  the  synop¬ 
tic  meteorologist,  i.e.,  the  forecast  behavior  of  typical 
synoptic  systems  and  the  forecasts'  comparison  with  conven¬ 
tional  technique.  It  is  hoped  that  such  a  summary  and 
analysis  will  permit  an  easier  evaluation  of  numerical 
forecasting  by  meteorologists  not  directly  engaged  in  this 
work;  in  particular,  with  the  aid  of  the  forecasts  oh  own 
in  the  Appendix,  it  is  hoped  that  a  "feeling"  for  th~ 
performance  of  the  relatively  simple  prediction  models  under 
a  wide  variety  of  synoptic  conditions  may  be  obtained. 

The  purpose  of  this  section  is  accordingly  to 
describe  the  data  preparation  and  to  present  the  synoptic 
study  and  summary  of  the  forecasts  made  using  the  simple 
barotropic  model  and  the  thermotrop.'c  model. *  The  synoptic 
study  of  the  forecast  results  includes  an  observation  of 
the  general  level  of  performance  of  the  numerical  techniques 
and  of  the  systematic  errors  which  are  in  evidence.  A  dis¬ 
cussion  of  the  forecasts  is  included  with  consideration  to 
initial  flow  patterns,  types  of  systems,  and  geographical 
areas  covered  by  the  forecast.  Case  studies  of  interesting 
situations  are  used  as  illustrations  of  performance.  Thc 
daily  forecast  and  observed  maps  for  the  500  mb  and  1000  mb 
levels  and  1000  -  500  mb  computed  vertical  motion  are  in¬ 
cluded  In  the  Appendix.  The  use  of  the  thickness  (1000  - 
500  mb)  for  forecasting  purposes  is  of  particular  interest 
to  the  synoptician,  since  there  have  been  ?  number  oi 
studies  made  in  the  past  to  empirically  relate  the  thick¬ 
ness  patterns  to  the  flow  patterns  at  the  surface  and  aloft 
for  forecasting  purposes.  Some  numerical  thickness  fore¬ 
casts  ar<  presented  here  in  detail  to  give  an  indication  of 
the  performance  of  this  portion  of  the  thermotropic  model. 
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A  discussion  of  the  thermotropic  model  as  a  forecast  tech¬ 
nique  is  included  and  a  comparison  is  made  to  other  methods 
of  prognosis. 

4.2  Synoptic  Data  and  Analyses 

The  decision  was  made  to  prepare  one  month  of 
carefully  analyzed  maps  at  all  the  standard  levels  to  serve 
as  a  test  series  for  the  numerical  prediction  techniques. 

This  would  provide  maps  for  sixty  sets  of  forecasts  which, 
it  was  felt,  would  be  a  large  enough  sample  to  determine  the 
level  of  performance  of  the  forecast  equations.  A  winter 
month  was  preferable  since  it  provides  the  most  difficult 
and  extreme  forecast  cases,  and  a  recent  year  was  preferable 
from  the  viewpoint  of  data  availability.  Examination  of  maps 
of  several  winter  months  resulted  in  the  choice  of  the  month 
of  January,  1953,  as  most  suitable  for  testing  purposes. 

Data  for  the  month  of  January,  1953,  were  obtained 
from  the  Daily  Upper  Air  Bulletin  published  by  the  United 
States  Navy,  the  data  files  of  the  Department  of  Meteorology, 
Massachusetts  Institute  of  Technology,  and  from  the  data  files 
of  the  Atmospheric  Analysis  Laboratory  of  the  Geophysics  Re¬ 
search  Directorate.  The  data  were  entered  on  WBAN-1  maps,  a 
Lambert  conformal  conic  projection  with  a  scale  of  1:12,500,000, 
for  the  following  eight  levels  for  0300  GMT  and  1500  GMT  of 
each  day;  surface,  1000,  850,  700,  500,  300,  200  and  100  mb. 
Figure  4.1  shows  the  areas  of  data  coverage,  analyses,  data 
tabulation  and  forecast  results. 

The  analyses  of  the  map  series  were  begun  with  the 
surface  maps  for  each  day,  with  time  continuity  established  by 
tracing  pressure  centers,  troughs  and  ridges  f  *om  map  to  map. 
Data  were  sparse  in  the  Pacific  Ocean  area,  so  the  coded 
analysis  of  the  San  Francisco  Weather  Bureau  office  was  used. 

The  upper-air  maps  were  analyzed  on  a  light-table  with  the 
analyzed  maps  for  the  next  lower  level  beneath  them,  which 
aided  the  maintainance  of  vertical  continuity,  although  dif¬ 
ferential  analysis  would  have  been  preferable  if  manpower 
and  time  had  allowed.  A  smooth  contour  analysis  depicting 
the  major  flow  patterns  was  desirable  for  numerical  predic¬ 
tion  computations,  and  this  was  easily  carried  out  on  the 
lower-level  maps  while  still  fitting  all  radiosonde  data. 

At  the  200  and  100  mb  levels,  however,  the  accumulative  temp¬ 
erature  errors  caused  sc  mingly  erratic  station  reports, 
which  indicated  very  irregularly  shaped  systems  and  seemingly 
unrealistic  flow  patterns.  The  analyses  at  these  levels  were 
therefore  in  general  tied  to  the  lower  layers  to  depict  smooth 
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large-scale  flow  patterns.  Frequent  errors  In  the  radio¬ 
sonde  data  received  necessitated  a  close  check  on  all  map 
data.  Reported  heights,  which  did  not  appear  to  be  con¬ 
sistent  or  were  not  in  agreement  with  the  general  flow 
patterns,  were  re-computed  from  the  ground  level  using  the 
equivalent  isothermal  layer  method.  Some  stations  consis¬ 
tently  reported  either  lower  or  higher  heights  than  sur¬ 
rounding  stations  due  to  the  use  of  different  radiosonde 
equipment  or  calibration  techniques,  and  this  was  taken 
into  consideration  in  the  final  analysis.  There  were  a 
total  of  496  maps  analyzed  for  the  month  of  January  1955. 

When  the  analyses  were  completed  a  square  grid 
of  17  x  23  points  was  placed  on  each  map  (Fig.  4.1).  At 
each  grid  point  the  height  was  interpolated  to  the  nearest 
ten  feet  and  entered  on  the  map.  The  values  were  checked 
and  transferred  to  tabulation  sheets  in  preparation  for  the 
computations  carried  out  on  a  high-speed  electronic  computer. 
Special  attention  was  given  to  the  data  from  the  1000  and 
500  mb  maps,  as  these  were  the  levels  used  to  test  the  per¬ 
formance  of  the  thermotropic  model  and  that  of  the  baro- 
tropic  model. 

4.3  Characteristics  of  January  1953 


Before  undertaking  the  synoptic  summary  of  the 
forecasts  themselves,  it  is  appropriate  to  discuss  the  gen¬ 
eral  characteristics  of  the  period  for  which  these  were  per¬ 
formed  .  Records  of  the  weather  of  January  1953  show  that  it 
was  not  a  "typical"  or  "normal"  month  from  the  point  of  view 


of  actual  weather,  in  spite  of  the  appearance  of  the  usual 


number  of  large  scale  disturbances 
tr ernes  in  the  Northern  Hemisphere. 
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It  was  a  month  of  ex- 
The  United  States 


experienced  one  of  the  warmest  Januarys  ever  recorded;  cold 
polar  outbreaks  were  usually  weak  and  of  short  duration,  the 
most  severe  cold  outbreak  occurring  on  the  15th  and  16th  in 
the  Great  Plains  and  Mississippi  Valley  area.  Precipitation 
was  twice  the  normal  in  the  Pacific  northwest,  southern 


Florida  and  New  England,  and  generally  above  normal  over  the 
rest  of  the  United  States.  The  strongest  development  of  the 
month  occurred  when  a  surface  low  formed  in  Texas  on  the  8th 


and  moved  to  the  southeastern  United  States,  accompanied  by 
heavy  rain.  This  storm  then  moved  northeastward  and  gave 
heavy  snow  and  ice  to  New  York  and  New  England  (see  Appendix). 
On  the  14th  a  cyclonic  development  on  a  cold  front  from  the 
Pacific  hit  Salt  Lake  City  with  heavy  snow. 


The  mean  500  mb  map  from  the  project's  analyses, 
when  compared  with  the  normal  January  500  mb  map,  shows  above 
normal  heights  over  the  western  United  States  and  over  the 
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north  central  and  northeastern  United  States,  with  below 
normal  heights  over  western  Canada  and  the  southeastern 
United  States  (Fig.  4.2).  The  major  trough  was  located 
east  of  its  normal  position.  The  zonal  flow  near  the  west 
coast  of  the  United  States  and  Canada  was  considerably 
stronger  than  normal  as  indicated  by  the  positioning  of  the 
centers  of  height  anomaly.  The  mean  jet  stream  location 
(200  mb)  was  along  the  northwestern  United  States  border, 
then  southeastward  through  the  southeastern  United  States 
and  then  northeastward  off  the  east  coast  of  the  United 
States,  passing  just  east  of  Newfoundland.  The  favorite 
surface  storm  paths  for  the  month  were  from  the  north 
Pacific  along  the  northern  United  States  border  and  from 
the  southern  Great  Plains  northeastward  over  the  Great 
Lakes.  Surface  pressures  were  below  normal  over  all  of 
the  United  States,  except  the  far  southwest. 

4.4  Synoptic  Study  of  Numerical  500  and  1000  mb 
Forecasts 

The  verification  of  prognostic  maps  obtained  using 
numerical  methods  or  other  techniques  is  a  difficult  problem. 
Statistical  methods  such  as  those  employed  in  Section  3  can 
give  an  estimate  of  the  accuracy  of  the  forecasts  but  many 
questions  remain  unresolved.  From  a  synoptician's  viewpoint, 
the  final  analysis  as  to  the  value  and  usefulness  of  numer¬ 
ical  methods  depends  on  the  accuracy  in  forecasting  displace¬ 
ment  of  pressure  centers,  troughs  and  ridges,  the  accuracy  of 
the  wind  forecast,  and  the  usefulr-  a  of  the  forecast  maps  in 
predicting  the  actual  weather. 

A  comparison  of  prognostic  maps  obtained  using 
numerical  methods  with  those  obtained  using  the  conventional 
field  techniques  is  one  cay  of  estimating  the  general  level 
of  performance  of  the  numerical  methods. 1°  Comparisons  of 
this  type  have  been  made  in  the  past.  In  1 95 1  a  set  of 
.enty-two  consecutive  24 -hour  forecasts  at  500  mb  from  the 
ixnearized  harotropic  theory  were  compared  with  synoptic 
forecasts  at  the  Atmospheric  Analysis  Laboratory  of  the 
Geophysics  Research  Directorate.  Correlation  coefficient  ■ 
of  the  forecast  versus  observed  changes  and  examination  of 
the  errors  indicated  that  the  two  types  of  forecasts  were 
of  approximately  equal  accuracy;  the  average  correlation  for 
numerical  methpds  was  +0.73  and  for  synoptic  methods  was 
+0.7-;.  Gates^'4  showed  a  comparison  of  a  numerical  forecast 
and  a  synoptic  forecast  at  500  mb  ‘which  indicated  that  the 
two  methods  were  of  comparable  accuracy  in  the  case  examined. 
The  sixty  numerical  forecasts  of  the  present  series  were 


compnv  ..I  '.'it!:  forecasts  obtained  using  other  techniques  and 
in  addition,  were  studied  to  determine  the  reasons  behind 
the  distribution  of  various  types  of  errors,  in  an  attempt 
to  discover  weaknesses  in  the  method. 

4.4.1  Distribution  of  errors 


Examination  of  the  January  1933  forecasts  revealed 
that  there  was  a  noticeable  geographical  distribution  of  the 
mean  height -change  errors .  Figure  4.5  shows  this  error  dis¬ 
tribution  for  the  too  mb  barot ropic ,  [300  mb  thermotropic, 
1000  mb  thermotropic,  and  the  1000  -  y00  mb  thermotropic 
thickness  forecasts.  These  error  maps  were  obtained  by 
subtracting  the  ae*  -al  TF -in\.  n  heimht  charge  at  each  grid 
point  from  the  fore  am  ---to.:'  n .•  F •'  *  change  for1  each  of 
i he  sixty  cases,  end  .  cry  •  F  "  'v.  atgr uraic  mean  error  for 
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errors.  In  western  Canada  the  tendency  to  forecast  5C0  mb 
heights  too  low,  and  thickness  values  too  high,  resulted  in 
1000  mb  height  forecasts  having  a  large  negative  mean  error 
in  that  region.  Deep  1000  mb  lows  were  erroneously  forecast 
to  move  into  western  Canada  from  the  Pacific  during  the  fore¬ 
cast  series. 

The  forecast  results  were  studied  to  determine 
whether  there  was  any  consistent  relationship  between  the 
forecast  error  and  the  initial  flow  patterns  at  500  mb. 

The  errors  were  plotted  at  the  grid  points  for  each  of  the 
sixty  cases  and  isoline  analyses  were  made.  The  geographical 
locations  of  centers  of  error  of  more  than  300  feet  were  noted 
relative  to  the  trough  and  ridge  pattern  at  500  mb.  Figure 
4.4  shows  that  there  is  a  tendency  for  positive  height  errors 
to  be  iocated  ahead  of  the  500  mb  troughs  and  the  negative 
errors  to  be  located  to  the  rear  of  the  500  mb  trough  and 
ahead  of  the  following  ridge.  This  is  evidence  that  the 
numerical  methods  generally  forecast  too  slow  a  movement  of 
the  trough  and  ridge  patterns.  The  same  impression  is  ob¬ 
tained  by  merely  looking  through  the  series  of  forecasts. 

Error  centers  adjacent  to  the  boundaries  were  not  used  in 
Fig.  4. 

The  numerical  forecasts  were  more  accurate  over 
the  eastern  United  States  than  elsewhere  over  the  grid.  This 
fact  was  shown  by  the  average  point  correlations,  root-mean- 
square  errors  and  average  errors  (Section  3).  The  day-by-day 
observed  height  changes  were  as  large  over  the  eastern  United 
States  as  over  the  rest  of  the  grid.  The  size  of  the  actual 
height  changes,  therefore,  cannot  be  used  to  explain  the 
difference  In  forecast  accuracy  over  the  grid  area.  Factors 
which  may  have  contributed  to  the  differences  in  accuracy  are 
the  hi *h  data  density  over  the  eastern  United  States,  the 
distance  of  the  eastern  United  States  from  western  and  north¬ 
ern  computational  boundary  errors  (which  are  the  boundaries 
most  affecting  the  forecasts),  the  inaccuracies  in  the  data 
over  the  -western  portion  of  the  grid  due  to  the  mountains 
and  the  flow  barrier,  and  thermal  block  effects  of  the  moun¬ 
tains  in  the  west.  The  high  level  of  accuracy  of  the  fore¬ 
casts  over  the  eastern  United  States  indicates  that  with  an 
increase  in  the  density  and  accuracy  of  upper  air  and  surface- 
data,  a  larger  initial  grid  and  Inclusion  of  terms  in  the 
equations  to  take  the  mountain  effects  into  account,  we  might 
be  able  to  Improve  the  numerical  forecasts  over  a  much  larger 
area  to  a  level  of  accuracy  approaching  that  attained  over 
he  eastern  United  States. 

The  500  mb  forecast  average  error  maps  show  negative 
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errors  in  the  northern  portions  of  the  grid  area  and  posi¬ 
tive  errors  in  the  southern  portions.  This  indicates  that 
both  the  barotropic  and  thermotropic  models  overforecast 
the  zonal  wind  on  the  average  during  the  month. 

*1.4.2  Forecasts  at  500  mb 

Messrs.  Hering  and  Mount  of  the  Synoptic  Section 
of  the  Atmospheric  Analysis  Laboratory  included  a  study  of 
results  of  the  numerical  forecasts  in  their  recent  evalua¬ 
tion  of  forecast  techniques.  They  measured  the  500  mb 
gradients  forecast  over  six  stations  in  the  United  States 
using  the  24-hcur  numerical  forecasts  on  fifteen  days  with 
identifiable  surface  cyclonic  centers.  The  average  wind 
vector  error  of  the  forecasts  was  26  knots.  They  note  that 
the  Air  Weather  Service  forecast  capabilities  study  found 
the  average  500  mb  wind  vector  error  of  a  synoptic  fore¬ 
caster  to  be  25  knots.  The  distribution  of  numerical  pre¬ 
diction  errors  is  found  in  Fig.  4.5.  The  forecast  area  was 
divided  into  four  sections,  ten  degrees  of  latitude  by 
twenty  degrees  of  longitude.  The  average  zonal  and  merid¬ 
ional  winds  were  obtained  in  each  section  by  averaging 
the  height  differences  on  opposite  sides  of  this  10  x  20 
degree  area.  For  comparative  purposes  the  errors  obtained 
using  persistence,  extrapolation  of  past  12-  and  24-hour 
change,  and  the  normal  January  map  as  a  height  gradient 
forecast  are  presented  along  with  WBAN  and  numerical  fore¬ 
casts  for  each  section  in  Tables  4.1  and  4.2.  The  wind 
component  error  of  the  numerical  24 -hour  500  mb  forecasts 
and  •-''f  WBAN's  56-hour  700  mb  forecasts  is,  on  the  average, 
about  equal  to  the  observed  change.  There  is,  on  the 
average,  a  noticeable  overforecast  of  the  gradients  by 
numerical  methods  for  all  the  10  x  20  sectors  as  well  as 
for  the  single  stations. 

The  displacement  error  of  the  500  mb  trough 
forecast  by  numerical  methods  was  measured  at  the  latitude 
of  the  surface  low  for  the  fifteen  cases  with  identifiable 
surface  lows.  The  average  error  of  the  24 -hour  trough 
forecast  was  l.t  degrees  latitude  (or  about  160  nautical 
miles)  for  fourteen  cases;  the  fifteenth  case  had  no  identi¬ 
fiable  trough  on  the  verification  map.  Vox'ticltles  in  the 
500  mb  troughs  for  the  fifteen  cases  with  surface  lows  were 
computed  U3ing  a  6*  latitude  grid.  The  24-hour  maximum 
vorticity  change  forecasts  obtained  using  numerical  methods 
were  compared  with  the  observed  changes  (Tsble  4.3).  The 
numerical  methods  forecast  the  correct  rign  to  the  develop¬ 
ment  in  every  case.  On  the  average  their  forecast  error  of 
the  vorticity  change  was  about  half  of  the  observed  change. 
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Table  4.3  --  Maximum  Vorticity  Measured  in  500  mb  Troughs 

Using  a  6-Degree  Latitude  Grid 


Run 

No. 

Verifies 

NWP 

24 -Hour 
Fcst . 

Observed 

24 -Hour 

Initial  Final 

24 -Hour 

Observed 

Change 

24 -Hour 
NWP  Fcst. 
Change 

NWP 

Change 

Error 

1 

2  Jan  55 

15Z 

42 

55 

30 

-25 

-13 

12 

2 

5/03Z 

28 

35 

25 

-10 

-  7 

3 

6 

5/03Z 

22 

48 

18 

-30 

-26 

4 

7 

5/15Z 

22 

20 

21 

1 

2 

1 

8 

6/03Z 

28 

18 

35 

17 

10 

-  7 

9 

6/15Z 

32 

21 

25 

4 

11 

7 

27 

15/15Z 

32 

25 

28 

3 

7 

4 

28 

16/03Z 

3^ 

40 

28 

-12 

-  5 

7 

45 

24/15Z 

30 

48 

38 

-10 

-18 

-  8 

46 

25/03Z 

38 

30 

50 

20 

8 

-12 

81 

27/15Z 

25 

25 

28 

3 

0 

-  3 

52 

28/03Z 

3^ 

35 

32 

-  3 

-  3 

0 

59 

31/15Z 

40 

38 

50 

12 

2 

-10 

6o 

1  Feb  53 

03Z 

48 

35 

50 

15 

13 

-  2 
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i  VECTOR  ERROR 
HR  ftOOMS  FORECAST 


This  indicates  that  numerical  forecasts  generally  underfore¬ 
cast  both  increases  and  decreases  of  the  large-scale  circula¬ 
tion  intensity  as  measured  in  such  a  manner.  The  results  of 
Hering  and  Mount's  study  of  the  50°  nib  numerical  forecasts 
give  added  indications  that  numerical  methods  are  presently 
as  accurate  as  other  short-range  forecast  techniques  at  50°  nib 

4.4.3  Forecasts  at  1000  mb 

Hering  and  Mount  in  their  forecast  evaluation 
studies  also  compared  the  1000  mbnumerical  forecasts  with 
forecasts  made  using  conventional  methods,  as  represented  by 
WBAN's  surface  prognostic  charts  for  the  month  of  January  1953 
Correlation  coefficients  were  computed  between  the  observed 
30-hour  surface  pressure  changes,  and  the  changes  predicted  by 
WBAN  for  the  January  1953  forecast  series.  A  uniformly  spaced 
grid  of  twenty  stations  covering  the  prediction  area  of  the 
numerical  forecasts  was  used  to  obtain  such  values  for  the 
WBAN  prognostic  charts.  Correlation  coefficients  for  the 
30-hour  ourface  forecasts  are  shown  together  with  the  coeffi¬ 
cients  obtained  for  the  numerical  24-hour  1000  mb  forecasts 
in  Fig .  4 .6. 

These  results  show  that  the  two  methods  are  of 
comparable  accuracy,  with  the  coefficients  for  the  con¬ 
ventional  forecasts  averaging  slightly  higher  in  spite  of 
the  longer  forecast  period.  The  root-mean-square  coefficient 
of  the  series  of  fifty -eight  WBAN  prognostic  charts  was  0.75, 
as  compared  to  0.o4  for  the  numerical  forecasts.  It  is  in¬ 
teresting  to  note  that  WBAN  and  NWP  curves  are  in  phase  with 
their  good  and  bad  forecasts  coinciding. 

The  1000  mb  forecasts  were  analyzed  to  determine 
the  accuracy  in  predicting  the  displacement  of  low-pressure 
centers.  Verification  was  initially  intended  to  include 
all  cyclonic  centers  with  identifiable  initial  and  sub¬ 
sequent  2-* -hour  positions  fall Ir.r  wlthiij  the  forecast  grid. 

Ho. .ever,  the  boundary  conditions  produced  a  marked  influence 
on  the  predicted  position  of  those  centers  located  near  the 
outer  limits  of  the  grid  area.  This  effect  on  the  evalua¬ 
tion  was  reduced  by  eliminating  all  eases  with  final  storm 
locations  ..lthin  two  grid  points  of  the  boundaries.  A 
total  of  fifteen  cases  were  available  from  the  January 
forecast  series.  The  observed  and  predicted  12-  and  2*- 
hour  positions  are  shown  In  Figure  4.7. 


N 

• 1 


lN3l9ldJ3oO  N0I1V13WK)0 
Fig.  4.6 


For  comparative  purposes,  a  similar  evaluation  was 
made  of  the  forecasts  prepared  by  WBAN  for  this  series  of 
fifteen  test  cases.  The  surface  prognostic  charts  transmit¬ 
ted  on  the  facsimile  network  during  January  1953  consisted 
of  a  30-hour  prediction  of  the  pressure  pattern,  plus  the 
forecasted  18-hour  positions  of  pressure  systems  and  frontal 
systems.  The  predicted  24-hour  position  of  a  low  pressure 
center  was  obtained  as  a  mid-point  of  a  line  connecting  the 
l8-  and  30-hour  positions.  These  conventional  forecasts 
cover  a  period  which  begins  three  hours  later  than  that  of 
the  numerical  forecasts. 

Comparative  results  are  shown  in  Table  4.4  and 
Figure  4.8,  The  observed  final  positions  and  tracks  have 
been  superimposed  in  Figure  4.8,  such  that  the  position 
error  is  given  by  the  distance  from  the  forecast  point  to 
the  center  point.  The  direction  error  is  defined  by  the 
deviation  from  the  horizontal  line,  and  the  speed  error  by 
the  differences  in  line  lengths  connecting  the  initial  po¬ 
sitions  with  the  center  point  and  initial  positions  with 
the  forecast ^positions .  The  overall  accuracy  for  this  test 
sample  is  essentially  equal  to  that  of  the  conventional 
methods  of  pressure  pattern  prognosis.  The  numerical  dis¬ 
placement  forecasts  show  a  striking  tendency  to  underestimate 
the  speed  of  the  1000  mb  cyclonic  centers.  The  predicted 
?4-haur  speed  was  too  slow  in  each  of  the  fifteen  forecast 
situations,  resulting  in  a  mean  algebraic  speed  error  of 
about  minus  three  degrees  of  latitude  per  day.  A  bias  of 
an  equal  amount  was  noticed  also  in  the  12-hour  displace¬ 
ment  forecasts.  This  systematic  error  appears  to  be 
primarily  due  to  a  tendency  to  overestimate  the  actual 
contour  height  at  500  mb  in  regions  of  anticyclonic  cur¬ 
vature  and  underestimate,  although  to  a  lesser  extent,  the 
contour  height  in  trough  areas.  Since  cyclonic  centers 
are  commonly  located  between  a  ridge  to  the  east  and  a 
trough  to  the  west  in  the  500  mb  pattern,  the  error  tends 
to  position  the  1000  mb  storm  center  too  near  the  trough 
or  cyclonic  vorticity  center  at  5°0  mb. 

A  comparative  evaluation  was  also  made  of  fore¬ 
casts  of  Intensity  changes  for  this  sample  of  fifteen 
cases.  The  measure  of  Intensity  was  arbitrarily  chosen 
as  the  gradient  of  pressure  over  a  six  degree  latitude 
distance  averaged  for  the  four  cardinal  directions  from 
the  center.  Results  of  the  test  are  summarized  In  Table 
4.5.  Figures  In  the  table  show  the  number  of  occurr  nces 
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Fig.  4.7.  itorm  tracks  (foracaat  and  obaarvad). 

Caaa  numbara  iacludad 


FORECASTED  CHANGE  IN  AVERAGE  GRADIENT  IN  MBS 


TABLE  4.5 

FREQUENCY  DISTRIBUTION  OF  INTENSITY 
CHANGES  OF  SELECTED  CYCLONES 


OBSERVED  CHANGE  IN  AVERAGE  GRADIENT 
IN  MBS 

(N)  W8AN  30  HR  SURFACE  FORECAST 
N  NWP  24  HR  1000  MB  FORECAST 
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Table  4 . 

4  -- 

Compa,  ison  of  Nur 
i  'recasts,  1. 

rical  and  Synoptic  24 
o  cases,  January  1955. 

-hr 

Direction  Er^or 
in  Deg. 

Spt  ed  Error 
in  MPH 

Position 
Error  in 
Deg.  Lat. 

Ave . 

Mean  Algebraic 

Ave.  Mean  Algebraic 

Numerical 
Forecast  - 

10  .7 

0 

8.2  -8.2 

5-4 

l\  BAN 

Forecast  - 

12 .1 

R5.2* 

5.8  2.0 

5.2 

*R  deslgna:es  deviation  to  right  when  facing  down.  ‘•roam. 
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of  a  particular  combination.  Numbers  in  parentheses  are 
the  results  of  the  WBAN  thirty-hour  forecasts;  other  num¬ 
bers  show  results  of  the  24-hour  predictions  by  numerical 
methods.  It  is  interesting  to  note  that  intensification 
of  surface  cyclonic  circulation  was  predicted  by  the 
numerical  method  in  only  one  case,  whereas  it  was  observed 
to  occur  in  eight  of  the  fifteen  situations. 

In  addition  to  the  study  of  the  surface  (1000 
mb)  forecasts  an  attempt  was  made  to  relate  the  accuracy 
(as  indicated  by  correlations)  of  the  numerical  forecasts 
to  the  initial  conditions  on  the  surface  map.  No  signif¬ 
icant  relationship  was  found.  With  initial  conditions 
which  appear  to  be  similar,  the  accuracy  of  the  numerical 
forecasts  varies  greatly. 

It  was  noted  that  in  the  immediate  area  of  low- 
level  cold  polar  outbreaks,  the  1000  mb  height  rises  were 
not  forecast  very  accurately  by  the  numerical  methods. 

This  may  be  due  to  the  fact  that  this  low-level  cold  air 
advection  is  generally  an  ageostrophlc  advection  and, 
therefore,  cannot  be  detected  by  the  equations.  The 
overall  correlations  do  not  reflect  this  error  as  the  area 
wherein  this  advection  takes  place  is  generally  small  com¬ 
pared  with  the  total  grid  area. 

During  the  major  cold  polar  outbreak  of  trie  month 
on  the  lpth  and  l6th,  the  averaged  correlations  of  forecast 
versus  observed  height  changes  for  the  four  cases  were 
above  the  monthly  average  at  both  500  mb  and  1000  mb.  The 
four  (}00  mb  cases  averaged  +  0.79,  as  compared  with  the 
monthly  average  of  +  0.74,  and  the  1000  mb  cases  averaged 
+  0.70,  as  compared  with  a  monthly  average  of  +  O.o.  .  The 
movement  of  the  major  storm  of  the  8th  and  9th  was  forecast 
well,  but  the  height  rises  connected  with  the  cold  air 
advection  to  the  northwest  of  the  storm  were  not  well  fore¬ 
cast.  The  case  of  14  January,  when  Salt  Lake  City  had  a 
heavy  snowfall,  was  a  case  where  the  r.umeri  "a  1  method  pro- 
due  d  a  poor  forecast  of  the  movement  of  the  low  center. 

Thf  center  was  forecast  to  move  too  far  north  and  the 
height  values  of  the  center  were  forecast  considerably  low¬ 
er  than  were  observed. 
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4.4.4  Location  of  jets  and  fronts 

How  well  one  can  forecast  the  location  and  strength 
of  the  jet  stream  (or  zone  of  maxlm-tst®  winds)  at  500  mb  using 
the  numerical  prognostic  maps,  is  undoubtedly  of  interest  to 
the  synoptician.  A  case-by-case  examination  of  the  500  mb 
forecast  and  observed  maps  (see  Appendix)  revealed  that  the 
maximum  gradient  at  500  mb  is  usually  well  placed  by  the 
numerical  forecasts  but,  on  the  average,  a  tighter  gradient 
is  forecast  than  is  observed. 

The  ability  to  locate  the  surface  frontal  systems 
on  the  numerical  1000  mb  prognostic  maps  Is  another  point  of 
interest.  Examination  of  the  1000  mb  forecast  and  observed 
maps  leads  to  the  conclusion  that  the  surface  frontal  systems 
generally  cannot  be  successfully  located  on  the  prognostic 
charts.  Attempts  to  locate  the  frontal  systems  by  placing 
them  in  the  1000  mb  trough  lines  results  in  gross  errors  in 
a  high  percentage  of  the  cases. 

4.4.5  Vertical  motions 

The  series  of  numerical  computations  has  provided 
computed  instantaneous  vertical  motions  in  the  1000  -  500 
mb  layer  for  the  initial  time  for  each  of  the  sixty  fore¬ 
cast  cases.  These  are  the  mean  vertical  motions  through 
the  whole  layer.  The  maps  showing  the  geographical  distri¬ 
bution  of  these  vertical  motions  are  included  In  the 
Appendix.  It  is  Interesting  to  note,  in  general,  the  com¬ 
putations  result  in  mean  ascending  motion  in  the  quadrant 
to  the  northeast  cf  surface  low  centers,  where  forecasters 
normally  expect  major  precipitation,  and  descending  motion 
to  the  southwest  of  surface  low  centers,  where  rapid 
clearing  is  to  be  expected.  This  places  the  centers  of 
ascending  motion  to  the  east  of  the  500  mb  trough  and  the 
center  of  descending  motion  immediately  to  the  rear  of  the 
500  mb  trough,  as  might  be  expected  .The  normal  rate  of  speed 
of  the  vertical  motions  at  the  centers  is  about  4  cm/sec, 
with  maximum  ascending  motion  speeds  of  12  cm/sec  during 
the  month,  and  maximum  descending  motion  speeds  of  8  cm/see. 

The  map  of  vertical  motions  for  1500  GMT  of  16 
January  (Figure  4.9)  show3  a  typical  distribution  of 
values.  Ascending  motion  was  occurring  northeast  of  the 
low  center  in  eastern  Canada  with  mr.ximum  values  of  6  cm/aeo. 
General  cloudiness  and  steady  precipitation  was  occurring 
In  this  area.  A  large  area  of  descending  motion  lay 
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southwest  of  the  surface  lov;  center  with  maximum  values  of 
4  cm/sec,  and  generally  fair  weather  was  present  over  the 
whole  area  of  descending  motion.  An  area  of  ascending 
motion  was  associated  with  a  surface  lov;  center  in  western 
Canada,  but  the  maximum  strength  of  the  motion  was  confined 
to  a  smaller  closed  area  of  2  cm/sec.  Variable  cloudiness 
and  scattered  precipitation  was  the  condition  in  the  north¬ 
western  United  States  and  southwestern  Canada.  In  the 
southwestern  United  States  there  was  an  area  of  descending 
motion  connected  with  a  surface  high  pressure  call  in  that 
locality  with  generally  clear  skies.  The  maximum  vertical 
speed  was  4  cm/sec. 

4.3  Case  Studies  of  the  300-1000  mb  Thickness 

Forecasts 


The  thermotropic  model  employs  the  1000  -  500  mb 
thickness  value  as  one  of  its  parameters  in  forecasting 
for  the  500  mb  and  1000  mb  levels.  The  1000  mb  forecasts, 
by  the  nature  of  the  thermotropic  model,  have  errors  in 
them  which  are  introduced  by  both  the  500  mb  equation  and 
the  1000  -  500  mb  thickness  equations.  Thus,  with  a 
perfect  thickness  forecast,  the  1000  mb  prognostic  map 
would  still  have  the  errors  which  ./ere  present  in  the  500  mb 
forecast.  Since  the  1000  mb  verification'  not  indicate 
the  degree  of  accuracy  of  the  thickness  forecasts,  for  the 
above  reason,  and  in  consideration  of  the  interest  of  the 
field  forecaster  in  the  use  of  thickness  maps,  a  few  cases 
were  chosen  in  order  to  make  a  closer  tudy  of  the  perform¬ 
ance  of  the  thickness  equation  and  to  attempt  to  discover 
possible  sources  of  error. 

Three  cases  were  chosen  for  study  where  the 
correlations  were  relatively  high  for  the  500  mb  forecasts 
and  quite  low  for  the  100U  mb  forecasts,  indicating  that 
there  was  a  major  error  in  the  thickness  forecasts.  A 
fourth  case  was  chosen  because  it  was  unique  in  the  series. 
In  this  case  the  500  mb  correlation  was  low  and  the  1000  mb 
correlation  was  relatively  high.  The  thickness  boundary 
errors  were  averaged  for  each  of  the  thirty-one  points  and 
plotted  on  a  map  covering  the  verification  area  (Figure 
4.10).  The  largest  average  errors  appear  l..  the  northeast 
corner  of  the  grid  area.  The  values  on  the  north  boundary 
averaged  1.5  to  2  times  the  values  on  the  south  boundary. 
This  probably  accounts  in  part  for  the  relatively  large 
errors  which  occurred  In  general  in  the  northern  portion 
of  the  grid  area . 
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The  dates  and  correlations  of  the  four  selected  cases 
are  listed  below: 


Initial  Time  R500  R1000 


Case 

1 

9 

January 

o 

o 

o 

r<~\ 

O 

0.62 

0.08 

Case 

2 

7 

January 

0300 

0.80 

0.32 

Case 

3 

13 

January 

0300 

0.76 

0.17 

Case 

4 

19 

January 

1500 

0.28 

0.84 

Case  1 

This  case  is  of  greatest  interest  in  the  north¬ 
western  portion  of  the  grid  where  very  large  errors  occurred 
in  the  thickness  forecast.  The  initial  500  mb  map  (Figure 
4.11)  shows  a  moderately  strong  trough  over  the  central 
United  States  which  developed  during  the  forecast  period 
into  a  deep  low  center  over  Alabama  with  generally  straight 
west-east  flow  over  the  rest  of  the  grid.  The  1000  mb  maps 
showed  a  low  center  moving  from  the  Mississippi  -  Alabama 
border  to  Georgia  during  the  period.  The  initial  map  had  a 
deep  low  off  the  northwestern  coast  of  the  United  States  and 
the  final  map  had  a  low  center  just  north  of  North  Dakota. 

The  predicted  1000  mb  map  agreed  with  the  ob¬ 
served  map  over  the  eastern  United  States  but  was  very  much 
in  error  in  the  northwest,  where  a  deep  low  was  forecast 
in  the  region  into  which  a  ridge  of  high  pressure  was  ob¬ 
served  to  move.  The  500  mb  predicted  map  had  a  contour 
pattern  similar  to  that  of  the  observed  map  but  the  gradient 
was  forecast  too  strong  in  the  northwest. 

The  500  mb  forecast  error  map  shows  that  the  over- 
Jy  strong  gradient  predicted  in  the  northwest  was  a  result 
of  forecasting  600  foot  excessive  height  rises  in  the  north 
central  portion  of  grid,  and  600  foot  excessive  falls  in 
the  northwest  corner  of  the  grid .  The  latter  error  may 
vfell  be  a  boundary  effect,  but  the  error  in  overforecasting 
the  anticyclonic  buildup  in  the  north  central  is  a  feature 
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which  appears  in  many  of  the  forecasts,  and  is  apparently 
inherent  in  the  equations  or  methods  of  solution. 

The  1000  mb  error  map  shows  heights  predicted 
1000  feet  too  low  in  western  Canada  and  800  feet  too  high 
over  the  eastern  Rockies.  The  thickness  error  map  shows 
the  amount  of  error  contributed  to  the  1000  mb  map  by  the 
thickness  forecast.  There  is  an  area  of  1000  foot  thick¬ 
ness  error  in  the  region  of  the  large  1000  mb  error.  A 
400  foot  error  near  the  northwestern  corner  of  the  grid 
is  of  the  same  size  as  the  boundary  error  and  is  likely 
a  direct  result  of  the  boundary  error.  The  boundary  error 
values  north  of  the  major  thickness  error  of  1000  feet  are 
smaller  than  the  monthly  average  values  of  boundary  error, 
and  cannot  very  well  be  considered  important  in  producing 
the  large  error.  Examination  of  the  1000  mb  and  850  mb 
charts  revealed  an  ageostrophic  flow  across  the  packed 
thickness  lines  producing  cold  air  advection  in  the  region 
of  large  error  whereas  the  contour  lines  did  not  indicate 
this  advection. 

The  initial  thickness  map  shows  a  pattern  very 
similar  to  the  initial  500  mb  map.  Prom  these  two  maps 
there  is  no  reason  immediately  apparent  for  expecting 
warming  in  western  Canada.  Antlcyclonic  thermal  vorticity 
values  along  the  southwestern  Canadian  coastline  due  to 
slight  antlcyclonic  curvature  and  wind  shear  may  have 
contributed  largely  to  the  forecast,  as  these  values  of 
vorticity  would  have  been  advected  by  the  equations  into 
the  region  of  large  error.  The  forecast  thickness  map 
shows  how  the  ridge  of  thickness  values  was  built  up  in 
central  Canada  and  the  gradient  was  decreased  in  western 
Canada.  The  observed  thickness  map  shows  that  there  was 
actually  a  decrease  in  thickness  acrocs  the  western  half 
of  Canada  and  a  maintenance  of  the  strong  west-east 
gradient . 


Speculating  on  the  possible  reasons  for  the 
large  thickness  errors  in  this  case,  it  appears  that  the 
following  factors  are  the  most  probable  reasons  for  the 
errors: 


(1)  The  nongeos trophic  flow  in  the  low  levels 
may  be  an  important  contribution  to  the  thickness  decreases 
over  the  period  by  advection  of  cold  air  at  these  levels. 
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(2)  The  forecast  equation  may  not  give  proper 
weighting  to  its  individual  terms,  since  the  coefficients 

of  the  thermal  advection  terms  in  both  equations  were  empir¬ 
ically  determined.  Theoretical  work  done  since  the  tests 
were  carried  out  indicates  that  the  coefficients  were 
probably  too  small. 

(3)  The  mountain  barrier  effects  may  be  impor¬ 
tant  in  this  case,  with  a  very  strong  west-east  component 
to  the  wind. 

Case  2 


This  case  is  of  major  interest  in  the  area  just 
south  of  Hudson  Bay  where  there  was  a  600  foot  error  in 
the  thickness  forecasts.  The  initial  500  mb  map  (Figure 
4.12)  shows  a  very  flat  trough  over  southern  Hudson  Bay, 
weak  to  moderate  west -east  flow  over  the  United  States 
and  a  low  center  in  western  Canada.  During  the  period 
there  was  little  change  in  the  pattern,  with  small  height 
rises  over  southern  Hudson  Buy.  The  1000  mb  maps  were 
much  the  same  at  the  beginning  and  end  of  the  period,  with 
a  deep  low  off  the  northwestern  United  States  coast,  and 
a  weak  low  over  the  south  central  United  States;  there  was 
an  increase  in  the  1000  mb  heights  over  Hudson  Bay. 

The  predicted  500  mb  map  has  a  strong  ridge 
built  up  through  the  central  United  States  and  Canada  and 
a  strong  gradient  across  the  northern  United  States  and 
soutnern  Canada.  The  correlation  coefficient  is  misleading 
in  this  case,  since  the  forecast  and  observed  changes  are 
pretty  well  in  phase  as  to  sign,  thus  giving  a  good  corre¬ 
lation,  whereas  the  equations  overforecast  the  size  of  the 
changes  by  a  considerable  amount.  This  tendency  to  over¬ 
forecast  height  rises  in  the  vicinity  of  ridges  is  Illus¬ 
trated  again  by  this  case. 

The  predicted  1000  mb  map  has  a  low  off  the 
north*.  -.  ..  coast  of  the  United  States  much  deeper  than  ob¬ 
served  and  strong  high  centers  southwest  of  Hudson  Bay  and 
over  Lake  Michigan,  The  1000  mb  error  map  shows  the  extent 
of  the  overforecast  of  high  pressure  southwest  of  Hudson 
Bay  in  the  same  region  of  maximum  500  mb  error. 

The  thickness  error  map  shows  that  sizable  thick¬ 
ness  decreases  were  forecast  in  the  Hudson  Bay  region  re¬ 
sulting  in  errors  in  excess  of  600  feet.  Boundary  errors 


were  very  small  to  the  north  of  this  major  error  and  cannot 
account  for  the  errors.  The  predicted  and  observed  thick¬ 
ness  change  maps  indicate  that  the  strong  warming  which 
occurred  in  southwestern  Canada  was  forecast,  but  show  that 
very  little  cooling  occurred  over  Hudson  Bay  where  strong 
cooling  was  predicted.  The  predicted  and  observed  thickness 
maps  show  how  the  forecast  of  strong  cooling  over  Hudson 
Bay  resulted  in  a  poor  thickness  forecast.  The  thermal 
gradient  along  the  United  States  -  Canada  border  was  fore¬ 
cast  much  too  strong,  and  a  trough  over  Hudson  Bay  was  fore¬ 
cast  to  be  much  too  deep. 

The  initial  thickness  and  500  mb  maps  indicate 
that  the  forecast  of  lower  thickness  values  over  Hudson  Bay 
was  probably  caused  by  the  advection  of  increased  cyclonic 
vorticlty  values  into  that  region,  and  not  from  advection 
of  cooler  air  in  the  layer.  The  initial  maps  indicated 
some  slight  warm  air  advection  south  of  Hudson  Bay  using 
the  500  mb  flow  and  the  thickness  lines  as  the  temperature 
field.  It  appears  that  the  advection  of  vorticlty  was 
weighted  too  strongly  by  the  equations  in  this  case. 

case  $ 


In  this  case  there  is  a  sizable  thickness  fore¬ 
cast  error  in  the  region  near  the  western  United  States  - 
Canada  border.  The  Initial  and  observed  500  mb  maps 
(Figure  4. 13)  show  a  relatively  unchanged  zonal  flow 
through  the  period,  while  the  surface  maps  show  a  buildup 
of  a  polar  high  in  western  Canada.  The  predicted  500  mb 
map  shows  good  agreement  with  the  observed  map,  while  the 
predicted  1000  mb  map  shows  a  deepening  trough  extending 
along  the  western  United  States  -  Canada  border. 

The  forecast  error  maps  emphasize  the  accuracy 
of  *he  500  mb  forecast  except  near  border  points,  while 
showing  a  600  foot  error  in  the  1000  mb  forecast  along  the 
United  States  -  Canada  border.  The  thickness  error  map 
shows  that  thicknesses  were  forecast  more  than  **00  feet 
too  large  in  the  area  of  maximum  1000  mb  error.  Other 
large  thickness  errors  are  located  near  regions  of  large 
boundary  errors.  The  observed  thickness  change  map  shows 
that,  there  was  cooling  in  the  layer  over  a  large  area  in 
northwestern  United  States  and  Canada,  wherens  there  were 
thickness  increases  forecast  in  some  portions  of  this 
region . 
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The  initial  and  observed  thickness  maps  show  that 
the  thermal  jet  moved  southward  from  western  Canada  into  the 
north  central  United  States,  whereas  the  forecast  did  not 
move  the  thermal  jet  southward.  The  500  mb  flow,  as  well  as 
the  low-level  flow,  had  a  component  across  the  thermal  gra¬ 
dient  so  as  to  cause  colder  air  to  be  aavected  southward 
into  the  north  central  United  States.  It  appears  that  the 
equation  weighted  the  advection  of  anticyclonic  thermal 
vorticity  into  the  north  central  United  States  to  the  point 
of  forecasting  thickness  increases  rather  than  decreases. 
Results  of  the  verification  studies  carried  out  on  the 
month's  data  indicate  that  the  forecast  equations  tend  to 
consistently  overweight  the  vorticity  advection  effects, 
and  underweight  the  effects  of  thermal  advection. 

Case  4 


This  case  will  not  be  discussed  in  detail.  This 
was  the  case  of  a  poor  500  mb  forecast  and  a  good  1000  mb 
forecast  (Figure  4.14).  The  500  mb  errors  were  small  and 
scattered.  The  low  correlation  indicates  that  the  forecast 
and  observed  changes  patterns  were  out  of  phase,  and  the 
major  error  apparently  was  one  of  not  moving  the  systems 
fast  enough  from  west  to  east .  The  initial  thickness 
pattern  was  almost  identical  with  the  500  mb  pattern  as 
the  surface  flow  was  weak.  The  thickness  forecast  was 
similar  to  the  500  mb  forecast  with  errors  of  similar  size 
and  location.  When  the  thickness  forecast  was  subtracted 
from  the  500  mb  forecast  the  errors  cancelled  and  the  re¬ 
sult  was  a  1000  ...^  prognostic  map  with  practically  no 
errors.  The  similarity  of  the  two  forecast  equations  is 
shown  here  by  the  similarity  of  the  thickness  and  500  mb 
forecasts  and  the  error  distributions  starting  with  sim¬ 
ilar  patterns. 

4.6  The  Thermotropic  Model  as  a  Forecast  Technique 

The  thermotropic  model  can  be  regarded  au  a 
practical  operational  forecasting  technique  considering 
both  accuracy  and  time  involved.  There  are  several  ad¬ 
vantages  in  numerical  methods  over  the  conventional 
field  techniques.  Numerical  weather  prediction  is  an 
objective  technique,  not  subject  to  errors  of  oversight  or 
the  physical  or  mental  conditions  of  a  forecaster.  The 
assumptions  made  in  developing  the  equations  are  at  least 
Internally  consistent.  While  their  solution  requires  no 
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knowledge  of  preceding  weather  in  forecasting  for  a  limited 
area,  boundary  conditions  are  required  which  may  include 
information  of  past  weather.  Numerical  methods  are  con¬ 
sistent  in  that  they  will  always  give  the  same  forecast 
with  the  same  initial  and  boundary  conditions,  and  offer 
the  mean  vertical  motions  in  the  1000  -  500  mb  layer  as 
a  by-product  of  the  forecasts.  Finally,  a  number  of  factors 
contributing  to  pressure  changes  can  be  more  exactly 
weighted  individually  and  collectively  by  numerical  methods 
than  by  subjective  techniques;  vorticity  advection  effects, 
for  example,  are  objectively  determined.  An  important  as¬ 
pect  of  numerical  methods  is  that  they  are  still  in  the 
development  stage  and  can  be  expected  to  show  continued 
improvement  with  time. 

Some  disadvantages  are  also  evident  in  this  and 
other  numerical  models.  First,  one  may  mention  that  numer¬ 
ical  methods  are  sensitive  to  errors  in  analysis}  the 
analyses  of  initial  data  must  be  smooth  and  accurate  to 
obtain  best  forecast  results.  Sparse  data  over  a  region 
are  more  of  a  handicap  when  using  numerical  methods  than 
when  using  che  conventional  field  techniques;  data  are 
generally  required  over  a  much  larger  area  than  the  area 
of  forecast  results.  Second,  arbitrarily  assumed  boundary 
values  produce  errors,  and  errors  are  Introduced  by  finite- 
difference  methods  of  solution  and  other  mathematical 
approximations.  Most  models  now  in  use  do  not  yet  have 
provisions  for  taking  the  terrain  effects  and  thermal 
source  effects  into  account.  Therefore,  additional  errors 
are  Introduced  into  the  forecasts  from  these  sources.  As 
far  as  the  equations  are  concerned,  pressui  systems  act 
the  same  over  oceans  and  continents,  over  warm  sources 
and  cold  sources.  Strong  cold  polar  outbreaks  into  the 
United  States  in  winter  normally  occur  with  a  strong 
age^strophic  component  to  the  winds  in  the  levels  below 
500  mb;  quasi -geostrophlc  models,  such  as  the  thermotropic 
model,  cannot  predict  the  changes  occurring  due  to  this 
cold  air  advection.  There  are  some  cases  in  the  January 
1953  series  where  the  numerical  forecasts  missed  strong 
surface  pressure  buildups  which  occurred  with  strong  low- 
level  cold  air  advection.  The  area  of  the  grid  which 
was  most  affected  by  this  type  of  error  was  Just  east  of 
the  Rocky  Mountains  along  the  United  States  -  Canada 
border. 
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4.7  Summary 

The  tests  carried  out  on  the  thermotropic  and 
barotropic  models  indicate  that  numerical  prediction  has 
developed  to  a  point  where  usable  forecasts  can  be  made 
on  a  routine  daily  basis  at  upper  levels  and  near  the 
earth's  surface.  Comparisons  between  numerical  methods 
and  conventional  field  techniques  using  the  same  data 
show  that  the  numerical  methods  are  consistently  as 
accurate  a  forecasting  technique  as  other  available 
methods.  The  thickness  equation  in  the  thermotropic 
model  apparently  forecasts  as  accurately  as  the  $00  mb 
equation,  but  the  overall  accuracy  of  the  1000  mb  fore¬ 
casts  is  somewhat  lower  than  the  500  mb  forecasts  due 
to  the  accumulation  of  errors  from  both  the  500  mb  and 
thickness  forecasts. 

The  barotropic  and  thermotropic  500  mb  fore¬ 
casts  Tere  so  similar,  as  shown  by  the  statistical 
study  of  Section  3,  it  must  be  assumed  that  (l)  the 
equations  themselves  are,  in  effect,  nearly  identical, 

(2)  the  methods  of  computation  Influenced  the  results 
to  such  an  extent  that  the  results  were  forced  to  be 
similar,  or  (3)  that  January  1953  was  a  month  in  which 
the  atmosphere  was  equivalently  barotrooic  throughout. 
Cases  such  as  those  of  8  and  9  January  make  the  third 
assumption  appear  unlikely.  Sections  of  the  maps  on 
these  dates  showed  very  strong  thermal  advection  and 
yet  the  barotropic  and  thermotropic  forecasts  were 
similar  in  these  regions. 

There  are  a  number  of  causes  for  errors  in  the 
numerical  methods  which  must  be  eliminated  in  order  to 
improve  farther  the  technique. 

(l)  The  mountains  in  the  western  United 
States  and  Canada  are  Important  as  a  flow  barrier  and 
as  a  thermal  barrier  between  the  cold  polar  source  in 
Canada  and  the  relatively  warm  Pacific  Ocean.  The  effects 
of  the  terrain  should  be  included  in  numerical  methods  to 
improve  the  forecasts  near  mountainous  areas. 

(?)  The  method  of  boundary  specification 
should  be  improved  in  order  to  reduce  ths  errors  intro¬ 
duced  by  the  computations  involving  incorrect  boundaries. 
The  results  of  these  tests  have  shown  that  the  boundary 
errors  are  occasionally  an  important  influence  in  the 
forecast  even  at  interior  grid  points. 
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(3)  The  quasigeostrophic  character  of  the  equa¬ 
tions  should  be  altered  in  order  to  account  for  the  effects 
of  strong  thermal  advection  due  to  a  nongeostrophic  win] 
component.  There  were  several  cases  of  polar  outbreaks 
occurring  during  the  month  which  were  misforecast  by  the 
model . 


(4)  The  method  of  computing  and  advecting  vor- 
ticities  apparently  produces  errors  under  certain  condi¬ 
tions.  With  strong  shear  conditions,  for  example,  there 
is  a  tendency  to  produce  too  great  a  height  change  down¬ 
stream,  which  is  particularly  noticeable  in  several  cases 
where  strong  anticyclogenesis  was  predicted  in  a  region 
downstream  from  a  region  of  strong  anticyc Ionic  shear. 
This  anticyclogenesis  was  not  observed.  In  some  cases, 
when  the  vorticity  advection  and  thermal  advection  were 
working  to  produce  height  changes  of  opposite  sign,  the 
vorticity  advection  was  emphasized  by  the  equations  with 
the  result  that  height  changes  of  the  wrong  sign  were 
predicted . 


(5)  A  predicted  1000  mb  chart  may  not  best  be 
obtained  by  combining  thickness  forecasts  and  300  mb 
forecasts.  Errors  are  accumulative  in  most  cases,  re¬ 
sulting  in  a  decrease  in  the  accuracy  of  the  1000  mb 
forecasts,  as  compared  with  the  individual  500  mb  and 
thickness  forecasts. 


The  fact  that  the  thermotropic  model  employs 
thickness  forecast0-  Is  synoptically  interesting.  Thick¬ 
ness  patterns  relative  to  storm  tracks  and  the  development 
of  systems  have  been  the  subject  of  many  synoptic  studies, 
including  a  comprehensive  study  by  Sutcliffe  and 
^orsdyke t and  interesting  applications  by  George and 
Owens.  ‘  The  numerical  prediction  project  has  used  the 
thickness  forecast  only  as  a  means  i obtain  a  1000  mb 
forecast,  and  the  verification  program  was  concentrated 
on  the  500  mb  and  1000  mb  forecast  .  The  only  thickness 
forecasts  evaluated  thus  far  are  included  in  the  present 
study.  Considering  the  general  interest  of  meteorolo¬ 
gists,  particularly  synopticians,  in  the  thickness  map, 
a  more  complete  3tudy  of  the  numerical  thickness  fore¬ 
cast  i 3  in  order. 
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The  continued  development  of  numerical  methods 
promises  to  produce  better  forecasts  in  the  future.  The 
mountain  effects  can  and  should  be  incorporated  into  the 
equations.  The  boundaries  can  be  moved  far  from  the 
region  of  forecast  interest  to  reduce  their  effects  or 
other  methods  of  solution  can  be  employed.  Continued 
research  should  provide  answers  as  to  the  best  methods 
of  vorticity  computation  and  advection,  the  inclusion  of 
nongeostrophic  wind  components,  and  the  best  levels  at 
which  to  employ  the  numerical  methods.  It  is  hoped  that 
numerical  methods  will  eventually  be  able  to  supply  the 
synoptic  meteorologist  with  an  accurate,  objective 
prognosis  c-f  the  flow  patterns  from  which  he  can  predict 
the  precipitation  and  cloud  condition  by  associated 
studies,  and  by  consideration  of  local  influences  on 
the  overall  weather  pattern. 

3  •  General  Summary  and  Conclm  ‘ on 

The  series  of  numerical  forecasts  described  in 
this  report  constitutes  perhaps  the  first  extensive 
study  of  the  barotropic  and  simple  baroclinic  nonlinear 
atmospheric  models  over  a  wide  range  of  synoptic  con¬ 
ditions,  and  lay  the  groundwork  for  the  further  system¬ 
atic  study  of  the  behavior  of  numerical  prediction 
models.  The  further  research  suggested  by  this  work, 
which  has  been  discussed  in  sections  3.6  and  4.7  in 
connection  with  the  statistical  and  synoptic  forecast 
evaluations,  stems  from  a  number  of  general  character¬ 
istics  of  the  results  and  which  may  be  summarized  as 
fol lows . 

(l)  Perhaps  the  most  striking  characteristic 
of  the  numerical  forecasts  is  the  high  degi'ee  of  similar¬ 
ity  between  the  barotropic  and  thermotropic  predictions 
for  the  pOCt  mb  level.  In  -  aoh  case  of  the  series  of  60 
forecasts,  made  from  the  initial  data  observed  every  12 
hours  at  03 Z  and  1*3Z  during  the  entire  month  of  January, 
1953*  the  barotropic  solutions  have  been  found  to  account 
for  the  position  and  orientation  of  the  large-scale  dis¬ 
turbances.  In  some  cases  important  differences  in  the 
forecast  amplitude  changes  of  the  disturbances  were  found 
between  the  two  models,  but  were  on  the  average  small  and 
not  significantly  different  on  the  statistical  basis  of 
the  correlation  between  predicted  and  observed  changes. 
This  evidence  reflects  the  quasi -barotropy  of  the  atmos¬ 
phere  on  an  average  basis,  and  strongly  suggests  that  the 
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baroclinic  component  of  the  flow  introduced  by  the  thermo¬ 
tropic  model  at  500  mb  is  in  gener  ..  relatively  small,  and 
that  for  24 -hour  forecasts,  at  least,  the  atmosphere  my 
be  considered  either  equivalent-barotropic  or  thermotropic 
for  the  purpose  of  predicting  the  flow  in  the  middle  tropo¬ 
sphere  . 


(2)  A  second  important  characteristic  of  the 
numerical  forecasts  is  the  marked  influence  of  the  lateral 
boundary  conditions  imposed  during  the  relaxation  process. 
These  conditions  not  only  seriously  altered  the  computed 
changes  in  the  immediate  vicinity  of  the  boundaries,  an 
effect  which  was,  of  course,  anticipated,  but  in  addition 
frequently  extended  their  effect  well  into  the  interior 

of  the  forecast  area.  For  those  cases  in  which  there  'was 
a  large  change  in  contour  height  imposed  on  the  boundary, 
the  predictions  are  in  general  noticeably  poor;  this 
circumstance  frequently  occurred  over  the  northwestern 
U.  S.  and  southwestern  Canada.  These  boundary  effects  in 
general  appear  to  have  affected  the  barotropic  and  thermo¬ 
tropic  forecasts  in  much  the  same  manner,  and  in  the  regions 
near  the  boundaries  probably  account  for  the  similarity  of 
the  forecasts. 

(3)  A  third  important  feature  of  the  numerical 
forecasts  is  the  high  level  of  both  statistical  and  synop¬ 
tic  performance  over  the  eastern  half  of  the  United  States. 
In  this  area,  the  24 -hour  forecasts  for  both  the  500  and 
1000  mb  levels  correlate  with  the  observed  changes  at  0.8 

-  O.p,  and  consisted ly  display  the  correct  position  and 
orientation  of  the  large-scale  disturbances,  as  may  be  seen 
from  an  inspection  of  the  maps  it:  the  Appendix.  The  fields 
of  computed  vertical  motion  also  bear  out  many  familiar 
synoptic  relations  between  the  vertical  velocity  and  the 
synoptic  pattern.  A  tendency  is  evident,  however,  for  the 
models  to  slightly  underforecast  the  observed  intensity 
changes  and  the  speed  of  propagation  of  well  developed 
dl sturbances .  On  an  average  basis  the  models  appear  to 
produce  better  forecasts  in  the  eastern  and  southwestern 
U.  S.,  with  poorer  average  performance  in  the  central 
plains  States.  This  average  behavior  is  suggestive  of  an 
important  effect  of  the  Rocky  Mountains  upon  the  flow,  and 
in  several  specific  cases  the  terrain-free  equations  were 
unable  to  predict  adequately  cyclogenesis  in  the  lee  of  the 
western  mountains. 
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In  studying  the  errors  of  numerical  forecasts, 
not  only  must  the  physical  approximations  introduced  in 
the  models  be  kept  In  mind,  but  the  errors  introduced  by 
the  numerical  methods  of  representation  and  solution 
themselves  must  be  considered.  Through  the  use  of  a 
specially  prepared  set  of  synoptic  charts  embodying  all 
available  data,  it  is  felt  that  the  purely  observational 
and  analysis  errors  have  not  seriously  affected  the  over- 
all  character  of  the  present  series  of  forecasts.  I;...  t tv- 
regions  most  likely  free  of  the  destructive  effects  of 
erroneous  boundary  conditions  and  of  the  effects  of 
mountainous  terrain,  namely  the  eastern  United  States, 
the  forecasts'  errors  in  many  cases  suggest  an  important 
effect  of  the  relative  vorticity  in  comparison  to  the 
coriolis  parameter.  This  effect,  it  will  be  recalled, 
was  omitted  from  the  thermotropic  equations  and  may 
partially  ax'v  in  the  observed  tendency  of  the  equations 
to  over-de;-..  lop  anticyclonic  disturbances  and  to  under¬ 
develop  cyclonic  disturbances;  in  the  atmosphere,  there 
is  observed,  of  course,  an  asymetrical  development  of 
these  disturbances,  .■.■hereby  the  cyclone  experiences  the 
greater  amplitude  change.  This  effect,  together  with 
the  systematic  truncation  errors  introduced  by  the  use 
of  spatial  finite-differences,  may  account  for  the  under- 
forecasts  of  the  speed  of  many  disturbances.  n  some 
cases  there  appears  to  be  an  appreciable  eri\  ,  especially 
in  the  1000  mb  forecasts,  due  to  an  inadequate  representa¬ 
tion  of  low-level  temperature  gradients  by  the  500-ao  ..b 

thickness  pattern;  these  errors  may  also  be  due  to  the 
neglect  of  the  so-called  "twisting  terms  and  to  the 
neglect  of  the  effects  of  the  vertical  advection  of  vor¬ 
ticity,  which  will  be  most  serious  In  regions  of  strong 
vertical  motion. 

From  the  presen:  tests  of  the  barotropic  and 
thermotropic  models,  one  cannot,  of  course,  expect  a 
complete  documentation  of  all  of  the  vorticity  generation 
and  transport  effects  operating  in  the  atmosphere.  It  is 
felt  that  the  present  tests  have  shown,  however,  that 
even  these  relatively  simple  models  are  capable  of  account¬ 
ing  for  a  large  portion  of  the  observed  behavior  of  the 
large-scale  atmospheric  disturbances,  and  that  in  areas 
free  of  the  more  obvious  sources  of  error,  perform  at  a 
level  comparable  with  if  not  slightly  in  excess  of  that 
achieved  by  the  conventional  synoptic  forecast  techniques. 


Much  further  research  remains  to  be  performed  in  order  to 
gain  a  clearer  understanding  of  the  behavior  of  large- 
scale  atmospheric  disturbances ;  research  is  in  progress 
at  the  present  time  to  investigate  more  fully  the  effects 
of  the  imposed  boundary  conditions,  the  effects  of  moun¬ 
tainous  terrain,  and  the  effects  of  the  use  of  finite - 
differences,  as  the  first  stages  of  a  program  to  isolate 
and  study  the  specific  effects  of  the  various  physical 
processes  themselves.  At  their  present  stage  of  devel¬ 
opment,  however,  the  present  techniques  of  numerical 
prediction  are  felt  to  be  suitable  for  operational  ap¬ 
plication  from  the  standpoints  of  feasibility,  economy, 
and  overall  reliability. 
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APPENDIX 


On  the  following  pages  are  presented  the  24-hour 
numerical  forecasts  prepared  from  the  thermotropic  and 
ba^otropic  models  from  each  IpOO  G.C.T.  synoptic  chart  for 
1  January  through  31  January  1953,  a  total  of  30  cases. 
Each  day's  forecasts  are  presented  in  a  series  of  eight 
maps  on  two  pages  arranged  as  follows. 

On  the  first  page  of  each  forecast: 

upper  left  -  the  observed  initial  500  mb 
map  (initial  conditions)  . 

upper  right  -  the  observed  final  500  mb 
map  (verification). 

lower  left  -  the  24 -hour  barotropic  500 
mb  forecast . 

lower  right  -  the  24 -hour  thermotropic 
500  mb  forecast . 

On  the  second  page  of  each  forecast: 

upper  left  -  the  observed  initial  1000  mb 
map  (initial  conditions). 

upper  right  -  the  observed  final  1000  mb 
map  (verification). 

lower  left  -  the  vertical  velocity  computed 
from  the  thermotropic  model  at 
the  time  of  the  initial  map. 

lower  right  -  the  24 -hour  thermotropic  100C 
mb  forecast . 


107 


>v  ->v»* 


S  .  >  « 
*\  ,  • 


l*  ,• 

\ . 

.  V  ' 


5«'K>  ve  ez  jan  <v-.« 


..  J 


l _ 


500  MB  AZ  5  J*«  l*>3 


f  X 


r*'<4 


k>o  mb  ©z 


8  JAN  1993 


500  MB  0Z 


9  JAN  1953 


'  .  «  "  -o  .  *  H  ‘ 


fcA»O?N0WlC  »OAICAft  »0  MS  ‘V  7  AN  *V! 


*  "*  *W0  1  ruAfCAtt  900  N*  »5Z  *  AM  (til 


1000  IO  152 


8  JAN  l<*5J 


WO  Me  02  to  JAN  1933 


300  MB  SZ  1  I  JAN  1953 


£xr  "\i\  c/«  7  n>  [ 

7,/  //V 

^  V  *y  \V\-'  •/ 


#7-:.  /■  /  y 

f  1  *;■  f>*  >y\\ 


f  ••  ■  -  *  ••  [,>  x  -V  •<•>*%•>  >■': 

.  .  .f  .  :  >.  V  ,  : « 

•AroTnonc  tom'-ift  non  at  uw  n 


»}/  m  .»  » 6  r  T* 


L. 


13  JAN  1953 


I0CO  MB  157 


f  ♦  .  ■  .^oe  .  c* 

i<X>0  M0  I3Z  14  JAN  1933 


1000  MB  IV  13  JAN  1933 


^vVn 

^  e  ..-  ? ,  -  j 

■  /  '  < 
/-  i  .‘'1  . J 

.  v  -1  .... 

tm  «f)  wj  ooo  Mt 


rl  j.  /  _  -  ^  >.  i  i 

vViXX'X  :^.v  A.'v  ■  M 

■  //%%!? ■< ; 

r\  yrr^ - -  •  ♦/•>•.>•<  xv- 

X-  Vi.-'.'-rr  -."V./.  V  .. 

'*  I  \  ,  ■'  /-■  ' 

•K  \ . ,••  • ;  •  • 

j  -V  •  ■  i  ■  / 

»  -  .  .  .  .0*  .  *  .  Of 

r*'*'c  coo  «)  »i 


ill 

m 

Vs 

0 

1000  MB  132  I  a  JAN  1999 


t  'x  ■ 

V'  .*  *  ;  * 


. /•  -  /°° 

O'  f  J  • 

^  ■  /  /  / 

•  •  -/■/.  / 
t  /  y-r-' 

'  V  •  v 

.  .  V  .N-'-rrT 

"**  '■**  x  - 

-  -■  K-X 

V'Ov  ....  V 


VVV, 

" m 


w 

/ A 

Ai  >  ‘ 

<  v- 


>■*  ••  .;»x  y  ;  x  t/" 

■  -v  •  •  •  f  ?  *;•  *  '•■.■  •  •-  ■ 

-  •  *  -  V  ■  •  .  ;  X 

' . ^s»<  ■>  .  .  q<L 

T«*B0tlWC  IIKtC *»»  U00  «  «  :t  M  *» 


19  JAN  1953 


1000  MB  152 


iooo  mb 


20  JAN  1953 


>■  v**  «*»rW 


500  08  BZ 


500  M6  ez 


2  2  JAN  1953 


t _ 


25  JAN  1953 


*•  «-  «* 


romcAj?  900  IM  t? ;  ;»  «) 


t*. 


.C 


*0 


,  %  V  \ 


l 


a*MtM0Y*m 


>  -•  )  / 

■  i 

'  / 

•  \  1  « 

.  \0f 

fOAtCAST 


r-  / 

I  * 

.  u  ;  4»  \  «c 
w  ip  *v 


*y  i  \  *•*> 


500  ue  02 


23  JAN  I9JJ 


•J 


•  •• 

500  “»  «  24  JAN 


'S3  3 


t-lftfctHQPC 


V  ♦ 


’  *■  .  * 


\#4 

W  \\> 

V  O-  v" 

'  t  V  v  ‘  v  ' 

■  ' 

,\* \  • 

*•  . 

'ONtcas?  V30  mi 


■  -  1 


f4  4441  0*> 


j 

\  • 

H-** 


w-^vvV  \  \  V  \  ; 


\ '  v 
.  '  \ 


;«  .  v*-  "\  ■.  .  \ 

k  \  \  \\  \  ■ 

w  v  y  \  \  \  \ 


\V  v 

kX  \  \ 


:>■ 


y.*t 


>\* :  :  '\n  V\\ \*\  V  • 

Ll1'\  XX\v\xn. 

\  W  x 


ft. 


■  V 


\ 


X:  X 

'  ■  ‘V  \  ■ 


vv>  ,  ' 

f  /  >  *>. 


\\\YV  •'  '//ft ■' • 


*  ■  v  \\\\,.*,'.  / 

■  •**.  ,\  .  ifc,  v  v  ■  •’ 


•  *4 

.4  MM  TO  J 


il*TiCA;  *C-T*C*  KM  UC 


500  MB  ez 


C3  JA*  1*3 


900  mb  ez 


\  V\  »o\  ** ; 


\  tq\ 

V\\\\  f 

V 

.v 


W  A\v 


\ 

x\' 


r  - 

*  .  ,  *  *  {  \  \  '  *■—  -•* 

•  t •  S  ••**'•  r  ‘  v*  ■  • 


. i 


(  * 
r  'y 


•■yv.ixv’-r :: 


\ 

’  .  K  K 

{  T*«WIOt»QnC  1 0*1 CAM 


wo  m 


'  ^<1  ffi* 


J  JAM  iWJ 


GEOPHYSICAL  RESEARCH  PAPERS 


Na.  1.  Isotropic  and  Non-Isotropic  Turbulence  in  the  Atmospheric  Surface  Layer,  Heinz  Lettau,  Geo¬ 
physics  Research  Directorate,  December  1949. 

No.  2.  Effective  Radiation  Temperatures  of  the  Ozonosphere  over  New  Mexico,  Adel,  Geophysics 
R-D,  December  1949. 

r>’o.  3.  Diffraction  Effects  in  the  Propagation  of  Compressional  Waves  in  the  Atmosphere,  Norman  A. 
Haskell,  Geophysics  Research  Directorate,  March  1950. 

No.  4.  Evaluation  of  Results  of  Joint  Air  Force-Weather  Bureau  Cloud  Seeding  Trials  Conducted 
During  Winter  and  Spring  1949,  Charles  E.  Anderson,  Geophysics  Research  Directorate, 
May  1950. 

No.  5.  Investigation  of  Stratosphere  Winds  and  Temperatures  Prom  Acoustical  Propagation  Studies, 
Albert  P.  Crary,  Geophysics  Research  Directorate,  June  1950. 

No.  6.  Air-Coupled  Flexural  Waves  in  Floating  Ice,  F.  Press,  M.  Ewing,  A.  P.  Crary,  S.  Katz,  and  J. 
Oliver,  Geophysics  Research  Directorate,  November  1950. 

No.  7.  Proceedings  of  the  Conference  on  Ionospheric  Research  (June  1949),  edited  by  Bradford  B. 
Underhill  and  Ralph  J.  Donaldson,  Jr.,  Geophysics  Research  Directorate,  December  1950. 

No.  8.  Proceedinga  of  the  Colloquium  on  Mesospheric  Physics,  edited  by  N.  C.  Geraon,  Geophysics 
Research  Directorate,  July  1951. 

No.  9.  The  Dispersion  of  Surface  Waves  on  Multi-layered  Media,  Norman  A.  Haskell,  Geophysics 
Research  Directorate,  August  1951. 

No.  10.  The  Measurement  of  Stratospheric  Density  Distribution  with  the  Searchlight  Technique,  L. 
Elterman,  Geophysica  Research  Directorate,  December  1951. 

No.  11.  Proceeding*  of  the  Conference  on  Ionospheric  Physics  (July  1950)  Part  A,  edited  by  N.  C. 
Geraon  and  Ralph  J.  Donaldson,  Jr.,  Geophysics  Research  Directorate,  April  1962. 

No.  12.  Proceedings  of  the  Conference  on  Ionospheric  Physics  (July  1950)  Part  B,  edited  by  Ludwig 
Katz  and  N.  C.  Geraon,  Geophysica  Research  Directorate,  April  1962. 

No.  13.  Proceeding*  of  ilia  Colloquium  on  Microwave  Meteorology,  Aerosols  and  Cloud  Physics,  edited 
by  Ralph  J.  Donaldson,  Jr.,  Geophysics  Research  Directorate,  'day  1952- 

No.  14.  Atmospheric  Flow  Patterns  and  Their  Repreaealation  by  Spherical-Surface  Harmonic*,  U.  Haur- 
witi  and  Richard  A.  Craig,  Geophysica  Research  Directorate,  July  1962. 

No.  IS.  Back-Scatteriag  of  Electrosmgaetic  Waves  From  Spheres  and  Spherical  Shells.  A.  L.  Aden, 
Geophvsica  Research  Directorate,  July  1952. 

No.  16.  Notes  on  the  Theory  of  l^rge-Scale  Disturbances  is  Atmospheric  Flow  With  Applications  tr 
Numerics)  feather  Prediction,  Philip  Duncan  Thompson,  Major,  U.  S.  Air  Force.  Gaopkysics 
Research  Directorate,  jsly  1962. 


GEOPHYSICAL  RESEARCH  PAPERS  (Continued) 


No.  17.  The  Observed  Mean  Field  of  Motion  of  the  Atmosphere,  Yale  Mintz  and  Gordon  Dean,  Geophysics 
Research  Directorate,  August  1952. 

No.  18.  The  Distribution  of  Radiational  Temperature  Change  in  the  Northern  Hemisphere  During  March, 
Julius  London,  Geophysics  Research  Directorate,  December  1952. 

No.  19.  International  Symposium  on  Atmospheric  Turbulence  in  the  Boundary  Layer,  Massachusetts  Insti¬ 
tute  of  Technology,  4-8  June  1951,  edited  by  E.  W.  Hewson,  Geophysics  Research  Directorate, 
December  1952. 

No.  20.  On  the  Phenomenon  of  the  Colored  Sun,  Especially  the  “Blue”  Sun  of  September  19^0,  Rudolf 
Penndorf,  Geophysics  Research  Directorate,  April  1953. 

No.  21.  Absorption  Coefficients  of  Several  Atmospheric,  Gases,  K.  Watanabe,  Murray  Zelikoff  and  Edward 
C.  Y.  Inn,  Geophysics  Research  Directorate,  June  1953. 

No.  22.  Asymptotic  Approximation  for  the  Elastic  Normal  Modes  in  a  Stratified  Solid  Medium,  Norman  A. 
Haskeil,  Geophysics  Research  Directorate,  August  1953. 

No.  23.  Forecasting  Relationships  Between  Upper  Level  Flow  and  Surface  Meteorological  Processes, 
J.  J.  George,  R.  0.  Roche,  H.  B.  Visscher,  R.  J.  Shafer,  P.  W.  Funke,  ff.  H.  Biggers  and  R.  M. 
Whiting,  Geophysics  Research  Directorate,  August  195. . 

No.  24.  Contributions  to  the  Study  of  Planetary  Atmospheric  Circulations,  edited  by  Robert  M.  White, 
Geophysics  Research  Directorate,  November  1953. 

No.  25.  The  Vertical  Distribution  of  Mie  Particles  in  the  Troposphere,  R.  Penndorf,  Geophysics  Re¬ 
search  Directorate,  March  1954. 

No.  2ft.  Study  of  Atmospheric  Ions  in  a  Nonequilibrium  System,  C.  G.  Stergis,  Geophysics  Research 
Directorate,  April  1954. 

No.  27.  Investigation  of  Microbarometric  Oscillations  in  Eastern  Massachusetts,  K.  A.  I'  lauraud,  A.  11. 
Meara,  b.  A.  Crowley,  Jr.,  and  A.  P.  Crary,  Geophysics  Research  Directorate,  Muv  195-1. 

No.  28.  The  Rotation-Vibration  Spectra  of  Ammonia  in  the  6-  and  10-Micron  Regions,  R.  G.  Breene,  Jr., 
Cap!.,  I  SAK,  Geophysics  Research  Directorate,  June  1954. 

No.  29.  Seasonal  Treads  of  Temperature,  Density,  and  Pressure  in  the  Stratosphere  Obtained  With  the 
Searchlight  Probing  Technique,  Louis  Kllerman,  Jul*  1954. 

No.  30.  Proceedings  of  the  Conference  on  Auroral  Physics,  edited  by  N.  C.  Gers.in,  Geophysics  Re¬ 
search  Directorate,  JuN  1954. 

No.  31.  Fog  Modification  by  Cold-Water  Seeding,  Vernon  G.  Plank,  Geophysics  Research  Directorate, 
August  1954. 
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No.  32. 

No.  33. 

No.  34. 

No.  35. 

No.  36. 

No.  37. 

No.  38. 

No.  39. 

No.  40. 

No.  41. 

No.  42. 

No.  43. 

No.  44. 

No.  45. 


Adsorption  Studies  of  Heterogeneous  Phase  Transitions,  S.  J.  Birstein,  Geophysics  Research 
Directorate,  December  1954. 

The  Latitudinal  and  Seasonal  Variations  of  the  Absorption  of  Solar  Radiation  by  Ozone, 
J.  Pressman,  Geophysics  Research  Directorate,  December  1954. 

Synoptic  Analysis  of  Convection  in  a  Rotating  Cylinder,  D.  Fultz  and  J.  Corn,  Geophysics 
Research  Directorate,  January  1955. 

Balance  Requirements  of  the  General  Circulation,  V.  P.  Starr  and  R.  M.  White,  Geophysics 
Research  Directorate,  December  1954* 

The  Mean  Molecular  Weight  of  the  Upper  Atmosphere,  Warren  E.  Thompson  Geophysics  Re¬ 
search  Directorate,  May  1955. 

Proceedings  on  the  Conference  on  Interfacial  Phenomena  and  Nucleation. 

I,  Conference  or.  Nucleation. 

II.  Conference  on  Nucleation  and  Surface  Tension. 

III.  Conference  on  \dsorption. 

Edited  by  II,  Reiss,  Geophysics  Research  Directorate,  July  1955. 

The  Stability  of  a  Simple  Buroclinic  Flow  With  Horizontal  Shear,  Lcoo  S  Pocinki,  Geophysics 
Research  Directorate,  July  1955. 

The  Chemistry  and  Vertical  Distribution  of  the  Oxides  of  Nitrogen  in  the  Atmosphere,  L. 
Miller,  Geophysics  Research  Directorate,  April  1955. 

Near  Infrared  Transmission  Through  Synthetic  Atmospheres,  J.  N.  Howard,  Geophysics  Res- 
search  Directorate,  November  1955. 

The  Shift  and  Shape  of  Spectral  Lines,  R.  C.  Breene,  Geophysics  Research  Directorate, 
October  1955. 

Proceedings  on  the  Conference  on  Atmospheric  Electricity,  R.  Holzer,  W.  Smith,  Geophysics 
Research  Directorate,  December  1955. 

Methods  and  Results  of  Upper  Atmospheric  Research,  J.  Kaplan,  G.  Schilling,  H.  Kallman, 
Geophysics  Research  Directorate,  November  195o. 

Luminous  and  Spectral  Reflectance  as  Well  as  Colors  of  Natural  Objects,  R.  Penndorf,  Geo¬ 
physics  Research  Directorate,  February  1956. 

New  Tables  of  Mie  Scattering  Functions  for  Spherical  Particles,  R.  Pesndorf,  B.  Goldberg, 
Geophysics  Research  Directorate,  March  1956. 


